For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


Ex AIGRIS - 
UNINEASTCACLS 


DIOS OD 

The University of Alber 
Printing Department 
Edmonton, Alberta 

a ee a ee 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Anderson1970_1 


7 = Note 
to tony f Si 
oa + ST eo . 


© ee arn riz 


U avPumere’. 
“SUBMITTED fo THE FACULTY OF GRADUATS aeons: 

aN BARTIAL PULPILMENT OF THE REQUIAEMNTS mua CUE ROMER 
: |. of pocron’ oF rurLasormy 


THE UNIVERSITY OF ALBERTA 


GEOMAGNETIC DEPTH SOUNDING AND THE 


UPPER MANTOE IN THE WESTERN UNITED STATES 


by 
@ CHESTER W. ANDERSON ITI 
A THESIS 


SUBMITTED TO THE PACULIY OF GRADUATE SIUDIES 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF DOCTOR OF PHILOSOPHY 


DEPARTMENT OF PHYSLcs 


EDMONTON, ALBERTA 


SPhiNG 60 70 


a ~~. 
Wa hai’ OUT aid, LT ‘eae 
emia TAY iit, ewes SP ian aD 


7 


* 
or 
- RE 
| 
: “= 
oe > 
ane F . 
| A 
‘ aa : 
_ 7 
4 
i a [ 
° 
a . a) 


: | arate), MPa 


ates anes 


UNIVERSITY OF ALBERTA 


PACULTY OF GRADUATES STUDIES 


The undersigned certify that they have read, 
and recommend’ to the Faculty of Graduate Studies Ae 
accepuance wa thesis ‘entuviled- GHOMAGNETIC DEPTH SOUNDING 
AND THE UPPER MANTLE IN THE WESTERN UNITED STATES 
SubmsL bed by Chester W.  Anderson UL an! partial fulfasiment 


Of ne Sequscements WOruie Cerree olf DOCTOR ml iitosopny 


a ¥- 7 re r ne i”, — 7 | | 
on . - i 
aati wr so eve 7 se din 
.. ROGINRS a2ayUeNe Ie ohh — 
We eat QETEM ORY. vin tn 
| aes ee ae ey alee 
jadgrnt afar eee are § ch haat : a bead | 


lain ant bi Mis ads300 to @Rkb old OT See 


ore 


ABSTRACT 


In September 1967 two magnetic variation events, 
a substorm and storm, were recorded by a two-dimensional 
array of 42 three-component variometers between latitudes 
So°N and 43°N and longitudes 102°W and 216°, Fourier 
aNealvels Ot these events shows that a complex spectrum 
with many maxima typify the storm and the substorm has 
aecmoOun spectrum. (Vpper mantle conductivivy soric cure 
Can De seen Gualitarively an the originals variograms bur 
is far more sharply defined in maps of Fourier spectral 
Component: “amplitudes “and phases. Separation of the 
FOUrTeEr= Spectral cOmpOnenvs Into internal and externa) 
Darts 21 Ves IUriner derintiiom tO Thewconduccive suruc— 
tures and allows quantitative analysis to be undertaken. 
PEP VOCeO fone COlduGtaVLTy Pills oom emae Dri NON rea ler 
than 200 km. under the Southern Rocky Mountains between 
therGreat Piains=and the Colorado Plateau, whieh marks 
aulOw—-Conduchivaty “recicn wathin toes Cordi itera. FA 
Strong conductivity anomaly runs north-south along the 
Waserch fronvy through central Utah, end andicates the 
presence of an upwelling of highly conductive material 
at sere no greater than 120 km. along the edge of a step 
Sserucbure Which Drings the conductive manule to shallower 


depths under the Basin and Range province than under the 
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Colerado platcau. The geomagnetic deep sounding anomalies 
ere found sto pe Bn excellent agreement with existing 
heav-1Jow data, thus supporting s correlation of electrical 
SOnaue rye ye WEL bemperapure.. Cbaere 1s also (200d 
COrrelation With the available seismic velocity informe— 
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INTRODUCTION 


Geophysical observations of several kinds indicate 
that the upper mantle of the earth under North America 
is laverally anhomogeneous. Upper mantle seismic 
velocities of 8.0 km/sec or larger are characteristic 
of Ghe Yeastern Unived States and the Great. Plains) province, 
Whe svelecrties tdecrease ico values sot 7.9 km/sec. or 
lower west of the Rocky Mountains (Herrin and Taggart, 
1962). SA Ssimivarwpattern 1s Shown by travel soime anomalies 
Gf S€lsiic Waves sav Vertical incidence, |) AS time-verm” 
like analysis of P and S waves (Cleary and Hales, 1966; 
Doyle vand Hales, 196/) gives evidence of early arrivals 
at stations in the eastern United States and predominantly 
lavesarrivals inthe western Unived staces, Asv che 
dirterences Wwevween P and S#uravel vime=residuals fand the 
eravicy anomalies in tae céentralsand westerm Unived States 
Canmousbe explained by the Birch (1961) réelavion between 
velocity and density, Doyle and Hales (1967) suggested 
thetetemperevure anomalies in the upper wmantle, possibly 
accompanied by parvial melting, might "cause lower 
velocities in the western United States by lowering the 
Valueston thevelastic constants. Archambeaus, Pilinn and 
Lambert (1969) substantiate the existence of partial 


Neibingipvarivaine evidence of) the eoexistence: of low-@ 
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ZOnes. Wathywones, of ow a velocities in the western 
United. otates. 

Lise nberpreval lon; ley supported. by eee heat 
flow data (lee and Uyeda, 1965; Roy, Decker, Blackwell 
encpe ech woos eB lackwell.. 1069 oewnieh, hear. Glow ante. 
Southern Rockies and the Basin and Range province 
COlre laues Wabi low Upper manvile vyelocitues:; normal heat 
Riove Gon LOUNGd. More tie Greate Piains, and, has high upper 
mantles Velocities, 

peismic and heat flow data are not abundant enough 
as yeu UO, delineate the boundaries, of different Upper 
Wein le provinces: Wilh saceUracy. | HOWeVer, .vemperavure 
eromalves im the upper mantle will) cause anomalies an tine 
Clecuraca mICOnOUCTIVLtye as sf bucates ate upper mneancHe 
temperatures show semi-conductor properties (Tozer 1959). 

“GeomagcnetL ie variations measured av the suriace: of 
che Marthe are the vector Sum of external vaerzatvions 
Orleamatiane. In Cumrente army tie  1onosphe re and magnevo= 
Spaere and of internal variations arising from currents 
induecdeby the external part. in the conductive regions 
ine the upper few hundred kKilomevers: Of tae Barth. lateral 
imnomnocenetples am the electrical conductivity of the 
upper manule Will give wrise to anomalies in, The, fveld 
Megcurvedsatmcne Surracel. — Qualitative anverprevagiony or 
Svcheconductbive Structures’ can be made from the amplitudes 
of geomagnetic variations in both the time and frequency 


domains if Fourier analysis of the variations is under- 
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taken, Gspecially-of the vertical component. Quantitative 
interpretation of such conductive structures requires the 
pepalavton Of the internal and external parts of the field. 
The amplitude and phase relationships in the frequency 
domaine between the internal and "external parts Gan then be 
titted by the appropriate conductivity models throuen the 
employment of conformal mapping techniques (Schmucker, 1969). 

In the summer of 1967 42 variometers designed and 
built by Gough and Reitzel (1967) to be inexpensive and 
peli@ble were sev Guu by the University of Alberta and the 
peoutmwesu Center for Advanced Svudies in a study “or the 
Ceusual and upper Mantle conductive Structures Im the western 
United States. 

Variation anomalies found during the summer of 1967 
Have Dee Oescrvped., Witte MuUalipaGive INverprevauuom ay 
heptzel, Gouen, Proratm and Anderson (1970)— A following 
paper by Porath, Oldenburg and Gough (1970) presents the 
results of a three dimensional separation of a polar substorm 
of September 1, 1967, and a preliminary interpretation by 
SLmMpD Le —cConductivauy models sor the rélavion between Une jup— 
SLOrimceparaved Lrelds.. “his Thesis presents lexvens 10m OL 
Vie Gua livavive Lntrerprevavion of UNG Variavion enomalies 
Dreseited Oy Reatzel eb al. 29/07 “Also Dreseited are results 
OlLwne- cevaravion of tne complicated polar svorm Ticia” oT 
September 20-21, 1967 with conductive models determined from 
the relationship between the separated storm fields. 

The writer joined the project in early 1967 and helped 


Sesemble and test twenty-two Gough-Reitzel variometers, prior 
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to the 1967 summer field season. The writer was a member of 
a veam consisting of Dr. D.i. Gough and Mr. W. Newman which 
insvealled. twenty-one Variometers in two east-west lines across 
the western United States. Regular servicings were carried 
out by Mr. Newman, Mr. McCullough and the writer at seven day 
Mover Vads-.) Tans “Servicing ‘routine required some "20,000 miles 
Om drive by cach menber O1ethe team. 

The weiter was LTesponsi ble rer the initial date. reduc 

tom Por the subsvorm of September 1),, 1967." The results of 

this analysis were presented at re 1968 American Geophysical 
Union mectingeurn Washineton, D.Ol, and av. the 1900) Canadian 
MSeSOCwatLOnrOl PhyclCiets Mesting. av Caleary, Witva. “The writer 
Wes respons ple for ell the data meduet ion and analysis of tire 
storm of September 20-21, 1967. Also the writer prepared maps 
Ch Uns eparaved Spectral “anplitudes mand phases of the: thiee— 
COMmponenus. OLS born the SUbSvTOLM and svorm tvelds mentwoned above. 
ines emresul uss were presented to the 1969 A.G.U. meeting in 
Washinevon. Di. anita paper oye heatzel ) Gough) rorach= and wnderson. 

TO eorroporave the model, or the upper manvule wnder vuhe 
western United’ States proposed by Porath et al.” (1970) tne writer 
seperated buree Components Of the svorm tveld in une trequency 
dotlain at? periods T=89, T=150 and T=256 minutes. Semi-quantitative 
iMcerprevarvom.oL thes storm antvernel fields aiethese periods 
deess support. Poratn's model’ particularly under tie Wasatch Fault 
Zone in, Uvahs 

RMaewwricver hasta lsouwshown ainvighapver 5 of This thesis 
that this model of the upper mantle is consistent with seismic 


ane heat tlow results in the western United States, 
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OUTLINE OF THE THESIS 


Chapter 1 develops geomagnetic deep sounding 
theory pertaining to the determination of conductivity 
Structure: underea plane surface starting from Maxwell's 
equations, through, to Seer ean estimates of depth and 
eonductivity structures in the Earth. The governing 
edUait ones lor clectromagnetic, induevion inp a conducting 
Hhaltespace wmulta—layéered structure and wndulating 
conducting surfaces are presented. Two methods of deep 
Sounding, vie "Base Stavion! and “Array” methods, are 
else presented. [he concept lor applyiney taste VourLter 
analysis to a single event. at each station over an entire 
2cray ss Gevel Oped Wit Lhe empnasts eon ume Physacal 
Meaning Of two-dimensional maps representing the sparial 
Vat PON Om ea MANS be Te Veni, meme aeaud OM Ol stile sb O@uat 
filed neon wes nternal end externals partss is. dts cussed 
Onna. shovel basiswand over a plane suriace:. 

Chapter 2 covers previous geomagnetic deep sounding 
imyvestucatsons In the wesvern United svaves.. Ihe work of 
Schmucker (1964), Caner et al) (1967), Gough and Reitzel 
G96] ) and Switt (1967) 1s reviewed with the aim of 
obtaining correlations with the work of heivzel et al. 
(197 0New Poravh et al, (19/0) (and this’ thesis as vo the 


locavions of the major magnetic anomalies in the western 
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UGITed obaves: primarily tees Waseuch Mront and Basterront 
anomalies. 
Chapter 3 gives a description of the 1967 array, 
the recording equipment used and the field procedures 
used., ‘Thissis Tollowed by aidiscussion on the processing 
OD thermagneticy Variation records from analom To dieital 
Porm,  eesectaon on the special spectral analysis 
techniques along with the method of separation of the 
Mapnerie Varlatvions is: also included. 
im Chapter 4, the basic analyses for the stations 
in the western United States from which magnetograms 
were obtained are presented and correlated with the 
mec swOolmnetezZelwet mal. ( MOVOD tendm ora coma len Wy 0) 
Chapter S discusses Unerusie Of Selsmolocy and neay 
flow work as upper mantle investigation techniques compared 
to the "Array" method of geomagnetic deep sounding. Sugges-— 
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CHAPTER 1 
GEOMAGNETIC DEPTH SOUNDING THEORY 


Mas “enapter develops @eomagnevic depth sounding 
theory aseappived 20 the determination Of conductivity 
SCCUC ure Under a plane surface, starting from Maxwell's 
cquatoie Winouch to theoretical estamaves: or Cep tn and 
CONGUCTIVIity ©f Structures in the? Earth. The governing 
euviav.ons= tor electromagnetic induction im a conducting 
half-space (flat-earth model), multi-layered structure 
ang undulating conducting surtaces are presented. Two 
Metnods Or dept sounding, the “Base wovativon” mernod 
and ytive® “Array method are considered. “The concept of 
apo yim eies 7 hOUrLe Vane y sais VOra siicle evemmmat seach 
Sve lOneOverwonte chulre array, is Ceveloped With the empnasis 
One thew Oinysileea limecanins Of vwo-dimens1 onal maps represen— 
Pine wUne we Dau la LAVarraApLOnS (Of saa Masnevrc even. woe para— 
Pr OMmOtmrne oOva le ie lds 2nv© it6 internal ene! eo weiimoer 
bocts Lerdiscussed on a @ilobal basis (Rikivake, Yokoyama 
andmooue, B9>6)"Walhinel and "Price 1939) and over a plane 
surface (Siebert and Kertz 1957, Weaver 1963, Oldenburg 
HI69 je Knowledge of Ghe “external field makes possiple 
SuUdvVaon tne MOorpnology of Une source Of the magnetic 
event. Knowledge of the internal field makes possible 


Be ceurave sLoOcavDi1on On Lhe map and estimation of maximum 
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depths of regions of enhanced conductivity within the 


Barth. 


Jes nel Narurenor the peomagnevic depth sounding problem 


Theta nducing TLtéld Tseusually of velopal dimers. ons 
OUG much Aanverest. lies in local variations of the magne-— 
Puce Lela Over "a lew tens or kilometers. The sphericivy 
Of the Barth can be neglected in such studies and the 
Hertha treated as a2 plane-surtaced infinitely large con— 
ductor having some non-uniform distribution cl Conducuivacy. 

EVE ligne. Inducing miei Beane bemva ken, as) mitiacorm 
Over onwalea,, Clee strengonh, On func, induced ~ucrentcrin 
Uwave area Wouncl, completely CStermimned by The sinducing 
fel ands thew eround "COncducrinvaty an the ammediace 
neighborhood The streneth or tne: inducedscurrenvs 
eepeniseone lie Cistriputaonyor The Senuire sen Cucan cam dew ad 
and whe seometry andy conductivity of a conduct om over 
ae reecien. OL similar dimensions as thaveot the inducing 
fieid (Dyck and Garland 19609). The geomagnetic depth 
sounding proplem then is the devermination if the re- 
Gistripution of some extended System of anduced currents 
Oy Local dititerences in conductivity strucvure over thay 


surface. 


ie ee Leto karth model with unuiorm Conducuivivy 


Mocts.01 wine inducvion theory presented Gere as 
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Wes L.o Meri ceniO5 0. 062 )F 
Por ene “rlat-earth" or intinive Nalr—space model 


Maxwell"s equations reduce to: 


VxEK = =e Cupei ide) tine VYxE = oar TMs ew iie 
VxH = 0 conductor VXH = 4nd Concucero@ 
Wolk = 0) Cz<0)) Wier S16 ease) 
Ve! =a.0 Ng DY 18 
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Be= cClecerical freld ingensuty Cvolus,metver) 
B = magnetic field induction (webers/m?) 
DD = dielectric displacement (coulombs/m?) 


H = magnetic field intensity (amps/meter) 


Qi eC mGmever spew la Lin CepencdenG mune Lome. 

ive siletace, Of they Conducvor as represenved by che 
mae OV omtanme= andy 2 as vaken posivive Cdowiwerd.) "by taane 
Pie VCMMaeoOt Nea cwiel Nome Cuatsons anda ies leculic tied s— 
placement currenus DOun inside and ourside the conductor 


Uiewreculvuine Wave equations reduce US 
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: 0 
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induction Equation) 


Tt Can be Shown that since Gg = o(z2) only, all 


Currents will flow parallel oO Une suriace Ole Uiesconductor 
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and Eo = 0 everywhere in the conductor. From y.J = 
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Lévy the electric fheld be of the form 


Hay ,2,0) = 22,8) xy) Oe 


where (1.3) becomes 
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where k? is a constant of separation. The governing 


equations Ansilde the conductor vo be Solved then are 
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for, the horizontal field and 
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for uhe vertical tield. Bavecion (1.0) also Wolds tor 


20. stor 2.0 the 2 equation is 
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‘ere SOlupLons DO the vertical field. 
Above che half-space (z<0) the external fiecid 
Can we Mritlen as whe @radient of 2 scalar potenrial 


@, which satisfies Laplace's equation. H = -VQ where 


Q = {Aexp(-kz) + Bexp(kz)}P(x,y,t)er”® : (io) 


(Prrce.. 190 7). 

| Physically the term A exe(-—kz) is the contribu- 
CLOn TO schice DOGettLal “Cue, UOneleCUrlc —CUurrenis 1imtume 
POnGspnere, and B exp (kz) is) the ongribuLion On une 
ineuCeG Currents wiulinm the conducroer= to, Une field ourside 
UViewconavetors. “ne med lis “Of §B/ Aw is tune ratio 1 Gmome 
emp lieude sO sue Taduced mMmagmep le. Trela TOs tvia Geom uiie 
imducame fLeld. whole are(B/A) as the phase datference 
DeuWeecnecne two fields: 

Price (1967) computed the effect of the source 
field in the table below on a "flat-earth" model by use 
Of Van LoeUcCt Mons parame ver (py — 4nuow/k? where o is the 
COMdUCtTI witty Of Une halt—-space., 7/2 Une frequency "or 
the Angueing field end K== 27/7, where ) is the wave or 
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fit seems then that induction becomes important only 
Wen Gol.) 8 2100 represents nearly complete control 
Or the induced current Dy self-induetion. The phase 
Guitoremces are Seen v0 decrease as 6 incredsc=,emeall 
Ohase Giitfevrences inply vnat Sselfr—induction 1s the 
dominant mechanism, while large phase differences imply 
Uier wmv cesistance ie the GCommnany, mecwanicm. 
Combining (1.7) with VxE = - 2] = - uH and using 


Tne parnonme Gime "dependence, Wie “exbernal wield becomes 


Hes a= fe CR ees pe ie 


oe 
Pie wPOUnGary Oni UVOns are bay pies vangcenul ale COmpOoneny 
Ol Ha2nG Ne norma COMmMpOnNens Of Bb be COntanmuous=acrocs 
the inverface. Applyine them to the solucvion of (1.7) 


Side conmparane leo) tol. So) produces he coc ft cient. 


nee ela he aie 2 Hw Gee) 
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MAGNETIC FIELD LINES 


FIELD VECTOR 
FROM SOURCE 


LINE CURRENT 
CURRENT INTO PAPER 


a ee 
\ RESULTANT / eS VITTTTTT177 
aN 
SX IMAGE SOURCE 
CURRENT OUT 
OF PAPER 


RIEED VECTOR 
FROM IMAGE 


Comparison of the coefficients ef 211) and (1. 12) shows 
Uheatecre ncuced horizontal components add to phe in— 
ducing field and the induced vertical component oynomtas 
the Inducing vertical fleld.) From Table l.i it-ean De 
Seen Chavet 6 > >, 8 > Avso tnav am the Limi ping case 
Of Very nigh conductivity Che magnetic field for any 


finste, weveleneuh as 
ES et ee ona Cis) 


ihe vhorizonval w;omoonents are then doupled white vine 
VeruLrcalaiLe Ga -componenvts are ance l led. 

The above result (1.13) can be illustrated graphi- 
Call yVacy sie wUse) Ol 1Meres lt Vanvexpernaly sine currens 
source is taken above a “flat-earth" model (infinite 
half-space) it produces an image of itself the same 
dispance below Utne boundary with en Opposilve Girecrion 


Of eUrremu st OW! (Hie wa. ty) 


133 Plane-layered model with uniform conductivity 


ihe wspavlal Conmuguratvlon of Chev surbacer mleldzon 
tne Z=Oeplane can be expressed as 42 spectrum Ol spatial 
Waveleneuns With @ coefficient representing vine amplitude 
CErcach. ~bet ouk)) be the ratio B/A (see section 1.2) orn 
the internal to external tangential fields at the surface 
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TAvcarauLe Wiki yary from Zerortor an insulator to one 


BOP To perteet. Concuctor. wet 
HES S.lG) y/ (ler tole) Gigs): 


(ie povential of a Magnevic Flela at the woundary 
Dlanew Cie. Mati —-space tai Tedspya\ plane avers. of var eus 


COMCUCTAVAU Nes Withee vine Last layer iexvendime to sits y 


(ice 1 Cai be expressed inte ditierent. corm wan 


equation (1.0) (Schmucker, 1969) 
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+ I(k,t)exp[i(k.r + kz)]}dxdy (1.15) 
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Mier veitdehinsa des the aes layer"is 
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=ai 
where ee = -k df _/dz aad ce Satasi Les The wdiiius.on 
equacvron Ul. (a). 
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ler eat dc= sie + (Hnwo,,)? 


ee then is the propagation constant of the downward difius— 
ing field in the men layer. By analogy with (i.7/b) 2 eeneral 


SOLUcHOMebOm eure ator bile ave layer is: 


tee) = A (Kt )exp(-k_z)+B (k,t)exp(k z) 


Wich 
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ie ra Oe G_(z)/k (CLABES) 
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An e - eee e 
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To \devermmine ~vhe-~vertical function, che continuicy 
COnGition for vne magnetic field vector ac the inter— 


beace A Ee is again used to give: the relation 


8 om 62 ima 
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m+1 oe m+ 


With Z=Z 4. equation (i520 when extpresscor liye mile ous 
hyperbolve functions yields the basic recurrence formula 
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m+l ~ “m 
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For most magnetic events used in geomagnetic-depth 
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The model wsed most frequently in the initaal inter— 


Oretation of depth sounding results 2s the Uwo-layer model 


Slpier Wwklokd UlLtOorm COnNGUCTLVIUy Or MOn-UniLorm Conducu1— 


Vitye hisssectilon well treat the uniform conductivicy 
@ace and show that dépth sounding and magnetotellurics 


are theoretically comparable in this limit. 
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The model to be treated can beuadescribed 2s 
uni rorm plane substratum of COnNGUCTAVILY OF that is 
covered by a poorly conducting layer extending from 
Si UNtO 45 24). The Ancident riela Genetrates 
UREOuUgi the Cop Vayer wie Mittle attenuation, whide 
Lis OeOUiNiOr Penelravion Minho Gie SubSstravum: tea sma lt 
compared to its spatial wavelength 27/k. Use of the 


relavLons 

ye Haley See 62 = i) / oe 
COsubseuLube Into equavilon ()\29) -erves. Mwave 7] Was 
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ine -compkex—-valuced paramever 
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Devine ravLo H7H or “/H We messured over ra sun— 
stratum then it is possible to compute its depth and 


Conductivity UnilgQueiy; thus 


h = w [Im (E/H) - Re (E/H)] , 


Gis) 


= 2 Re (E/a) 


1S 
W 
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for the two-layer case. 

1p is seem wien thay une real part one (out—o1— 
phase #/H or Z/H) represents the mean depth of the inter- 
Malveddy currents, while the imaginary pert of «¢ (in-phase 
E/7 nor 4/H)) yields with I Che amb lenny COMmMaue Inv yea 


that depth. The height 
hel = shee Cis) 


Gane mcnoughy aot sas Une “Aepun Of a perrtecucUbsciruce 
Conoucitor » Mo multi“ wavyered substravum Ts sindistinguisn— 
aolemiron, nd Uiereicre replaceable Dy. 2) Unit ern 
COMCUe TOR au thay depen no aoe ees oo eres lone eho 
Sinelewirequency component of the incidens Variation 
field 15 concerned. This equivalence makes the two-layer 
model very useful in the interpretation of deep sounding 
results in the frequency domain. 

Stnce F/H and 2/7 musth always be POSluive or wero, 
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eonductor" at various frequencies for 


Lahiri and Price Model d of the Earth: 
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Pie soup eL—plase component of B/N or 2/H will @iye mone 


G@epth to a perfect substitute conductor 
% awl 
h =w Im (E/H) (Age 


while the in-phase component gives the apparent conduc-— 
tivity 


Oo. = biSn(Re Ey) ]2) (187363) 
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atu thav depth. Fig. (1.3) from Schmucker (1909), shows 
Be VEresUsmOsLOrra*variepy Or Trequenctes=torstae Mani 
encwr Mice model yd concueclivivy Glscri bution: 

The Macnevoreliiuric derinivton of whe “acpareny 
Pesistavity (Cagniard (1953)) is based om the modulus 
Or une Ampedance and would yield as the apparent. skin-— 


depen value; the modulus of cy. 
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This section will ideal wrth the eChiecu of. laveral 
Conductivity Inhomosencviies Upon the surface magnetic 
fieldevariablons.. [hss uype Of SUrUuCctUre Can De wdelineared 
Veryeweld by the G.DsS. method with mespecy to ite strike 
andwlateral extent. amd to some degree tO ive cepun. 

The following treatment 1s restricted) 10 jessentially 
a theoretical representation of the first order effects of 


lateraleconductivity inhemogeneities on theveurface Picid: 
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Hence only the two-dimensional case need be treated. 
Since the anomalies will be considered in two dimensions, 
only one component. of the horizontal field and the 
Wertieal Tieldayill be mecescary to  descripe one induc- 
tion process completely. For ease in going from theory 
VOyiMverprevatton Of resulis in) later ‘chapvere, Z will 

be taken as positive down and y is parallel to the change 
PVsCOnauecy Lvalby, "rad oO. 

Two distinct types of anomalies will be treated even 
Laoughy Other types Of anomalies and combinations of 
anomalies exist. The surface anomaly is due to super— 
ticielecongucuavity Varketions in the cruse. Swthese 
enone les are Or, tNverest Tipper many re ives ead ul Ome 
Magnily because they can shield lateral conducciyity 
changes im he Upper Manvle “and Cause Misinverprevavsous 
Of Upper Manele Conductivity Ssumucture., “Tits siielda ns 
effect seems to have been in evidence along Caner's (1967) 
profile which was operated to map the northward extension 
Of the hao Grande anomaly. “Im the Wanta Besin anwUvan 
deep conducting sediments cause phase anomalies without 
appreciably attenuating the effects of upper mantle 
Structures below (Reitzel, Gough, Poravh and Anderson 
(eT) Chapver 4 of this thesis). “Ihe NowmtaeGerman 
anomaly is now known to arise from shallow sedimentary 


conductive structures (Vozoff and Swift 1968). 
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Figure 1.4 
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The second anomaly to be treated is produced by 
lateral conductivity inhomogeneities in the upper mantle, 
which can be représented by a perfect or nearly perfect 
Conductor withean undulating surface. .This second type 
of anomaly 1s of direct importance to the upper mantle 
Lav SS Utes Or. 

Either type of anomaly following Schmucker (1969) 
Can be Unought: of as the resultant of two ditferens kinds 
Of tavoma tous ~ivtetds. ihe dirst Kand arises from one 
PAOUG EOI Ol VOCAL “Currents. and sproducy ron, Ob aoe Oc laved 
eliomeavousm ir ernal Tielas Dy Whe mainly horazonval 
“normal” ‘field at the structure, and 1s represented: by 


the transfer functions h (Schmucker 1969). The 


H? “H? 9H 
SeconGden iG esreautus Irom normal anducvien ia the hori zon— 
Telly wayered parts of the structure by the anducing Zz 


field and is represented by hy» 27s Age 


Surface Anomalies 


The model considered (fig. 1.4) has a surface sheet 
end ean underlyine cConducvive Nalt-space separated by an 
ineulevine Weyer. (ne “surrace Can be Prepresenvesd bya 


Hon=UnvGorm woinm sheet with 42 verlable Total conduc sivivy. 
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where Jes the integrated conductanee in a sheet of 
thickness d, 
a 


d =f *o( 2) az 
O 


Ts thes” > "O plane., “Prices (1967) spoundary cond ons 
for the anomalous part of the variation field above (+), 


and below (-) a thin sheet give 


ol +t, _ = - ee ee 5 eS 
dy CH ce He )-(H - H) = iecepan Aniw[JZ,+ J,2] 
(1237) 
tte CWOMOUMeNnS ONS) sac 7 Us the X-cComponent, of stne mornoily 


ImaurCceQmeheet cCurreny denslyy for ee and ane Ge Cenotred 


in normalized form 


SS aye Gee ei) 


Fquacion (1.39) can be replaced by the approximation 
Die = Nie a in.) 


where Tiles TS vies percinenly 1nduccion sparvame ver ibs taken 


in one, source tela tree laniteis 
ie Ntwoc 


and the complex-valued c is given by equation (1.31). 
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The superimposed anomalous part of the induced sheet— 


CURrentmdens tly /i6 


cf A = 
dows (CH - ES 


where Ns is the result of both the firsso anda second 


Kinds, of antvernal freid 


ee = dH + dye : Vis) 


Equation (1.39) can then be expressed as 


= [(Q + dy )/J] = eniwzy, 
Creo) 
s (a,/3) = 2wiw(z, + J,/3) 


Hora gine enoma lous variations of the first and second kind 
respectively. The boundary conditions, equation (i.40), 
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A relaxation method suggested by Price (1949) 
when applied to equation (1.41) will give the induced 
field above an area where the conductivity distribution 
Weer The “surrace is. welt known. Whe stratified deep 
COMGUCUIVELY Slructure enters the problem a6 a tree 
Darameter CO be adjusted GO yea Che best iv, Gor ome 
observed anomaly. Price suggested two apoproximacion 
methods tO Solve the problem; 1) to negllect selr— 
imduction, 2) tbo meglect the Ohmic resistance. These 
eppreximMavr ons. are. comp lemenvary, sOumlav if sunemiarse 
COcsnOg Lead co COnVeEreing approxima bvons.swhes second 


Was ad On so. 


Undulalvineourtaces 


The reaument presented is an the mains tha toL 
Sehmucker (1969) and Rikitake (1965). 

Meavera | cCOnduGuiIVivy AnhDOMmoOrenelvies im Une upper 
Mantle cause: a Vara plewdep tn ae een i Ome ace 
TO piace for a Given frequency of 42 maenelic Variation 
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tangential to itss surface 
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dy 


= (Z + ey a Hae CW) 


ihe Tami by of field Mines that satisfy the’ boundary 
condition (1.44) and do not intersect the earth's surface, 
Z-0, form the possible interfaces between non-conducting 
end. perieculy conducting material. The fiela dines muse 
also satisfy the condition that they merge with a pre- 
Conceived Mocmal level Tor Gne perreco subsotcube conduc— 
Or, ee Some cvsvance trom the anomaly. hn! can oe 
computed from the equations developed in sections (1.3) 
end, (la). 

Wie sdowmweard extension Of Ue Vari atl omer oul dp its 
Carraed oul separavely Tor ics normal “and snomalous) panus. 
Love anomalous Dar Really as Oo) sinvernal sora cil. emer) 
ies downward exvension 18 Straiehntrorward. Tne anomalous 
field fan be expressed av z=0, by a series of Spatial 


hermonuce. tiea, Simusoidal field dtsvrelbuvion, 
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Also any attenuation by a surface conductivity layer 
iS disregarded, that is the shielding effect upon H 
de smolts” The elope of the field dine Jay = ay) 
Passing through the point (yz) follows from equerion 
Os 

Ae ee = ao coy 


= 2 (eG) 
ei) + Q'H 
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Ge envire set or field lines ¢an be put togevher 
Dy verapircal or numerical mevpnods. sstarting wie wi Cy )= 
fh! Se coint cistant from the anomaly. The ciame Lunar tome 
of equation (1.46) must be roughly in phase, meaning that 
Che se Gvenvatcon ine new sUrhaces Layers 1s “siaeas regu liced 
iio llivy eine ereumenv om Ol, seinem iby Vann (one) 4s 
Chic sisi coupared Wao Ulnmoy.  Theretore ius srype sot 
Dverprevatlon applies UO Cases Where sell —-imductioen a3 
dominant, that is = deep conductors of Larse radius) and 
ich conductiv ivy = 

Sehmucker (19609) considers a sinusoidal dnduetion 


anomaly ot the form 


rol 


Ge 0) = C sin Cae 
Z(Vs ee —C cos(k.y) A Cee 


with the conduction dH/ay = 0 
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Rikitake's model (1965). 
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Substitucing equation (1.47) ante (1.46) gives 


CBCOS WMOey))) exon z) 
ee 2 ge Steno see ee (1.48) 
apy co xe sin(k, y)exp(k,z) 


rol 
I 


exp (k,2) 


The transformation v=exp(-k,Z) allows equation (1.48) 
LOMDe solved analytically... The solucion: 


oy) Su aT sin(k,y) exp (k,z) . (1.49) 
a 


iihete ie alalkadaiay k2<<l the sine term will predominatre 
PiVine sanusoidal osciliarvsvons., while the: Lamat k2>>1 
PrOUUeces OS CMM atITOns 2 vine morn Ol sqvere waves. slme 
Field Times @oscillate asymmetrically about 400). 

Rikitake (1965) used a different method to study 
laveral conductivity inhomogeneities from the standpoint 
Of Si musOlos ) Undmlavwvons mf a Derrecu Subst ruLT =] scondue— 
won. Hewacsumed, The sUriace Of Mie "conductor to have a 
sinusoidal dependence and then analysed the etrect of an 
micident tield on the unduleavions , whiereas  camucker 
assumed = Sinusoidal anomalous field and, Convinued Unis 
fTelaucounwarcc) GO LOormuleve he mavure of Vie nda laving 
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Rikitake (1965) let the surface of a eonductor 
Be fepresented my Utne function 2 = f(y) on whieh the 
POvenvIal Of the Anduced magnetic field can be defined 


eca(see Gre, * 5), 
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Mesassumocion Of asurhace fly) = N= cos y leads toa 
sey Of Simulvaneous Equations which Can be solved [cr 


C 
n 


M4 
> 
D 

>) 
= 
ae 
J 
(ep) 
Il 
ae 
2 
3 
= 
O 
Ly 
Ey 
Il 
as 


(ep): 


LOrm aN ae 


le~al 
5 
Ga 
2 
> 
a 
> 
~ 
Ht} 
S 


39 


_ 


gel’ i= f 


. ' - i im 

; ao 
nt on iw ‘Ss nent Pape 
~ betuand aid toe nh 


t@nl +a}? os 


_ 


" 
j i o¥ 
i 1] 
i 
i e - 
’ > 
‘ 
’ } ; ny 
, < 
: 
. 1s 


yy 
¢ « - 
a. Logiwets, SF 
_ 


pm a es Bit: 


HO 


where 
" rh cos 
SON nD ee =) Y sin ny san Ny ay 
=i 
Pe nnicos 
an = Of 26 Y sin Vy “eos ny Sane Ny ay 
-T 
Lie et ne COSit Telents eonverzse and lve e mield 


CLStrlvouUoLon tSamiular Go that Or ochmucker within tae 
prescribed Limits Of ratios ‘of the field variacrion to 

tne Meprent Ol che undulatuons; bh ="Olla and k2<<l 

are ‘equsvaleny Jintts thet cive the Same sinusoidal 

Oc AMAL SCVOlS . AK vake  s) resitlps. ter bDaewundu apaor 
problem give “a Maxamum or the horizonval anducing field 
apove. che undulation and “a slighvly asymmcurie snduced 
Verulcad tveld wich a transition frome poOsitLve to Metavive 


direciity “above the=undulataon. 


in ‘Sones lees 


Pemvurbal Los. Of Tne magnevosphere by sexvernal 
sources produce eurrents in the ionosphere and. magneto— 
Splere that cause changes in the Barth's magnevic) field. 
These Magnevic Variations can be characterized by the 
Vector P1e1d.,, 10s time durataon andiits geograpnic exvenu. 

Schmucker (1969) shows that polar magnetic sub= 
storms have frequency components and amplitudes such 


that the depth of penetration is limited to 600 km. 
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Currents above this Limiting dépth @ive magnetic field 


componenvs Which Gan be detected on the Barth's surrace, 
Making events of this type ideal for upper 
Mave Studies. Subsvorm Ticlds are of Varge extent and 
COnscquently “sulvable Tor study cf structures across «4 
Continent. The morphology of Substorms (Rostoker 1960) 
ie KNOW CO SUCH “al extent as CO Dermay one Go make 
LOUgi tte Oreulca LV esvimetes. Ob. Lie  CmMaraCrerm lor siiae 
exverne |) imelG Of fa particular vevent. 

Bonnevier et.al. (1969) suggest that the parts 
OP UMe CuUecreniys, AsSsSOClaved Wilh polar Magmeule supe vormts. 
Waosewt Pelds sare ImpOrtany au une Karull a surimace,, can 
DEwGeOresenved Dy ta Three -dimens onal cumreny mode 
(fae leo) The cUrrents are: “Fa currently tlowme 
Gowniward ailene a> geomagnevic field Vine inve vae Lono— 
Sphere. 2) an anvense weseward electrojeo, 3) a revurn 
fieltd-=alivened currenv. Ine closed loops exvending ouu— 
Ware from uhe westward Electro jel cam be vaken as 


parallel Vo ve Harti Ss surtace and roughly along a 


geomagnetic equipotential. The gradient across these 

equipotenuials represenvs tne external horizontal) maz— 

MeuLle Variavion Observed On Uhe Bervi"s surface. 
iiewiine Pape lledaiD = 0. tien lO, m ic ule 


demarcation line of a substorm (Rostoker, 1966). That 
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Figure 1.6 


Possible currents in a polar magnetic 


substorm. Brcer Oldenbure. (l909). 
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POSSIBLE CURRENTS ASSOCIATED 
WITH POLAR MAGNETIC SUBSTORM 
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is, stations east of this line will observe a negative 
D bay and “stations wést of this line will observe 4 
posiwive DV bay. —“Rostoker (1960) has shown that if the 
current system itself moves, some stations will observe 
a transitional bay, i.e. the demarcation line moves 
Prom Cast vo west across the stations. 

(He -curreny model predvevs thay voe characrer— 
MSUReS Obea Subs LOI Can be expected to vary Wann 
Ceomagneric® lavitcude, as has been observed for meny 
VetLew we new aurora Zone Wherenr Meweleocuroyeu sins 
dominant. Larvernieretive—-m Days Will soccur Wiki poslelve 
JCoOmoOnents, CoO. ches norvly and negative 4 vo Une souch or 
(ee eu OraNele Cur oy et. 

In mid—-latitudes positvive-H bays and negavive-—Z 
bays are vo jbe expected af the perturbatrons are caused 
only by ine auroral veleccurojeu. (bso. unewiivensuty 10% 
DhHe Supsvorm should decrease as we proceed scuth of the 
Clrreny sysvem. The horizontal field “should also become 
More Nol rortiw —Excepuplons vO this General suvuatton occur 
during. partacularily intense Storm activity, when te is 
possible for the current system to move southward and 
brozgen itself to such an extent that observations at 


mid-latitudes show a southward increase of H (fig.4.12). 


inva Storm southward inereases of H may arise from the 
iveld- of che Pring-—current. in. the maznetosphere...bubr in 
mid-latitudes this increase will amount to only about 
five gammas.. 

Rosvoker 1968) hast @tated that the onset wf a 
polar magnetic substorm is always accompanied by at 
leas0 CWO Olsevines Pie micropulsations which cecur 
almost Simulwaneously at stations hundreds of kilometers 
Qari Us ile. Onset. bimes eo within S0lseconds) shor 2 
substorm can be obtained from magnetograms which clearly 


show these micropulsatvions. 


ie( lwo basic geomagnetic depth sounding mevhods 


Tae wo Maimwmeuhiods Om GDS Miivestaee ut Onmare 
wine Wase suation and, the “Array metunoda.. | [ne whiten 
GOsUeon magne vomevers and [ield programs  limived most 
Of the early work in GDS to the “Base station" method. 
Schmucexer (1964) Jand Caner et 421, (Cl96/) carried out 
extensive programs in the western United States with 
Dio sete COC s 

The sdesiom and ions tructlon Of a. 1oN-Cos Gg vhree 
component magnetometer by Gough and Reitzel (1967) 
permusoced the “array” method to be employed in the 


western United States. THe Writer was 2 collaborator 
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fivieoti est ticldeprogram emp loyineernis method. 
Forty-two magnetometers were used in the western 

United states during the summer of 1967 (Meitzel, 
Genghy Porach and Andersonsl9/0)¢ The array method 
allows simulvaneous recording of “the same variation 
eVenUs Over an area, with concentrations of magnetometers 
Over suspected anomalous Zones. A further advantage of 
une Varray metniod 1s that 1 makes" possible a formal! 
Separacaon Of PuUNnBRee. Componenvls Olea Vertar lon: 1emd 
TMVeOulnvernalyand Vexternal@ parce soy une use Vor ssuriace 
integral methods developed by Vestine (1941), by Price 
and Wilkins (1963) and by Weaver (1963) and applied by 
Clidenbure (1969), Porath, Oldenbure and Gough (1970) 


and the writer. 


Base station method 


Geomagnetic depth sounding aves Ula VOM Oye vile 
lipase Isbauion Mevmod Drocecac Dy recording Ginterent 
Magmetiec events Of the same type, such as Days and) sub— 
SvOrMs,au Gitferent locations ima specified (area where 
a-permanent observatory is located or where a portable 
magnetometer remains fixed during several events. The 
analysis of the magnetograms is then directed toward a 


statistical correlation between the anomalous and normal 
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Remo fou Tmesetevents.. Wieder survey, e@tatvionde ox 
matrix of transfer funcuLons Ls. obtained, connecting 
the components of the anomalous and the normal magnetic 
Vaeietion Vectors at that site in the frequency domain. 
Sehmuckéer (1969) states that these transfer functions 
WE elve Whe proper Stat wswical basis form thescibce— 
CUeWwemimMverprevation Of Cie imduct fon anomaly, Im8 perms 
OF tnvernal sconduchiyiry, surucvures. 

Tne “spaculea lh coniseuravion onthe primary snicucing 
Prelate sessenvially wagnored by vis met node Wares 
thus limited to mid-latitude and equatorial areas where 
Ge wex cernal wteld jean De assumed Wo We moze vor plees 


DG mern ene My Poravntetral. ©lbO70) Wievenenonny gies venwin 


mMmid—Lagpipudes the exvernal field can be far from uniitorm. 


imen Base wuacwon meted 1s derinively Bamived wan: 

determining the anomalous and meormal fields, ssimce, vais 
enalysis isa Tuncewonsor vne devarvhed KnowledZe sok ine 
llapeesecale variation treld whichwmay have seale Length 


Speorcer, 43000Ckm. 


Meray ame vliod 


Whes“array” mevnod proceeds by locauine Many 
magnetometers in a specified area (in our case 42 


magnetometers in a 700 km. by 1,300 km. area) which 
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remain fixed at their sites for several magnetic dis-— 
CUpbanecea. ihe investigavorwnay allow his station 
Spacing LO vary as he wishes, but the location of. the 
Station with respect nye en SLCUCCUrPeS Camilo us le 
as erDitrary. Other geophysical and geoloesical data 

Cale DewuUsedstO CeLermine “Zones Oo: possible anomalous 
SUBSUPUCtIUre. Thereby allowing cConcentracions of SstauLons 
inevnese Zones VO Obtain as Much Aniormation as possible 
about the suspected anomaly. The ideal method is to 
Operate an exploratory array fLirsy vo devermine the 
1ocavlon Of the anomaly, then (To Wwse a more concentraved 
array in the anomalous zone in a second anveswleoa yon. 

iin eseperetlvon as Fo We avueipued: -COnur Ol sbavLons 
aQuelaree ev ances feOm me anomaly. prenerapl yn Reem 
areas, Should be Operated Go give 2 good) Trepresentavion 
Of the external field with respect “tor 10s; scaler tenevh.. 
Scale lengths of magnetic variations in mid-latvavudes 
are uyplcally of the order 7,000 km: (Sehmucker 1964, 
PCrarmmep eal, 9 O70). 

The first order result of the “array” metned as 
aPewo-dimensional pilevure or a Single magnetic destur— 
bancesover-on area Of large eExtvenvu. ~The disturbance can 
be Fourier analysed atmeach station to give its specvral 
components. The amplitudes and phases of spectral com- 


ponents when plotted on maps give a two-dimensional 
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Pplevure on Une Vartarion aoress the array av a 
pParvicular Trequency. Section Wo willl 26 Into more 
Gevamlapoul such Mourler tranetorm maps. Whew array” 
memnog "also shows that much additional informariron can 
be “obvained "from the anomalous horizontal t2elds immduced 
Oyeuhne Normal "exvernal horizontal fields over Local 
invernal conducrors, wlth ae SUriKe Of Une WCOnduGceor 
alone cle Ete tne Change or wradieny =o» wne anomalous 
BOLT ZOntTa Pei Velo. sboun tne Wormzontalwand verulea. 
components map boundaries of the lateral conductivity 
PiMOmMoOeener lL hes Wile 2reaver “accuracy when mosey mevlhods 


Gr Upper Manole wiavesulcau Lon. 


1.8 Fourier transform maps 


Lae din ws sone ef Une external Magnetic (ved imo 
Coe pearr CMe CWitUSnOn wor Foie e tue ldsmOn sCie sii he imal 
eddy euUrrenvs and’ tac induction process im a conductave 
Sseuructure are aid Gependent "on whe Brequency on Fue 
ex vernal. magnetic Variation field. ihe reliavmony becween 
tne external and incernal fields, whose vector sum forms 
PhewVvarIavilon iaeld Observed eal UNG: Suriece, is fa1so 4 
funecion of frequency. ~ Therefore more TihOrmtat Lone Can, De 
obtained from the amplitude and phase relations av par- 
ticular frequencies than can be extracted from the 


original date an the time-domain. Thus the concept, OL 
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a vwo-dimensional representation of these spectral 
Components at a@ particular frequency,  , can be 
CNOUgNT Obeas @ direct consequence of the “array” 
Concep. Doward beuter determinations of the Earth's 
eonductivaty sveructure. This use of Fourier transform 
maps Was developed by Reitzgel, Gough, Porath and 
Nndersom (19/0 )" 

pe limaved areavolr wne Hartin Ss Supmace sean pe 
Pepresenved Dy avplane, 2Z=05 Upon which observations 
OMVUne smarMetic Variavion fteld at Many Svarions are 
to be Made 7) Spectral estimates of aycomporent, <(t); 
Ome une ey atl ay vont Held ‘ean De sderermimedwal weachn om 
UNneSewSvava Ons a Accordance wir vie equations 
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On replacing the invegral wath a sul Of Tinive terme; 


N-1 


eta eer pa an cos (21k) + Dd sin(Hk)) 
J O i k N k N 
k=1 
z 
[eo ec) m2) 
Were aj =a0,l sc, 16. eNaL eave 2 real numbers 
Xy> Xa» X55 cae Xoy_y> Which. cOMmpuces. Une our lei re kine 


and Gosine series coefiricients Ags ay> Di> Ass» Dos 


a b See The relations 
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el oere, 
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eas 
Arcten i) = 6). 


give the amplitude and phase of the magnetic variation 
MG) sate el particular frequency component k, respectively. 

fOr a) Particularstréqueney ¢omponents kK yiune 
values Ay and OL ap each Station can ber plotred ena 
map Vand “Contoured Using linear wnterpolation as an 
approximation. 

if one. camension of “an anomaly is motezreaver 
thane the corresponding dimension of the array, then 
Chive anomaly wilt produce, Sradivents am the anplatude 
ao unabecani De anterpreved Mala tetivelyes iim werine so ft 
CONGUCtTAyUUy Strucoure. “These eradients, would be: due to 
inipernal| currenus: ravhier Than exvernal sourcesmvecauce 
the external sources vary fairly smoothly across a mid— 
latitude array (section 1.0). Although the ismoothness 
Oiarwinewmnelds Gecreases, Toward the auroral welectiro; cis, 
the gradients produced by major conductivity changes 
Mey solid be CGistanguashable insmaps of unmeeparated 
Fourier transform parameters. 

The phase’ is) complementary to the amplivude? in 
the same manner as the sine and cosine transiorms are 
complementary. Thus an anomalous field gradient may be 


observed in the phase maps while little or no anomalous 
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behaviour is detected in the amplitude maps. The phase 
maps can contain information that may confirm interpre- 
tations made from the amplitude maps. An example of 
such a confirmation follows. Lateral conductivity 
inhomogeneities can often be modelled by a step structure 
Ghateshould produce an asymmetry in the Z amplivude maps 
ecross ene step structure. If anJadditional hump is 
placed on the edge or Chis step then a reversal should 
OGCUMelT wae Zecomponent., Since she phase maps contain 
Information about the sign of the anomalous field, the 
Peversal Ou. should sappear as “Steen eradienve 2m ecnes4 
Phase wmaps.s lhws was "exactly, tne case whem tne Woscavcm 
fhrony anomaly was sbeine TniGerprered from emolirude and 
phase unseparated maps (Chapter 4). 

the samolitude and »phase unseperaved maps allow 
ayeocd times inverpreteavilon ol anonalous Welds ii werms 
ObmCOnGUCtian CY WeLrucvure Cla qualitative lasier 

Phe numerical magnitudes of whe sine and icosine 
transforms (sine and cosine Fourier coefficients) have 
Chysical meaning only in relation to the sign and 
megnibude: of bhe spectral component Concerned av sche 
Starts oteche: time Inverval, ©, transformed § —Achange 
imeune jscarbing time will change the pellative: magnavudes 


and signs of the cosine and sine transforms, and the 
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GiUsetrPoOuUrionm Ol 4 e@iven spectral term of the time 

series between the cosine and “sine coefficients’. 
TmGrcder VO eseparate the Total field aAnto ats 

external and internal parts in the frequency domain, 

tie seine ane cosine transform maps of each must be 

separaved, Independently bDecause the Surface integrals 

Usecw tor (seperevicon can only be applied) to fields in 

eonstant phase across the map. The separated components 


com tertaver recombined) Dy une re lav. ons 
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where ae and PF, are the external and anvernal amplivudes 
and o. and On are the external and invernal, phases, 
respecvively., 1 The Variavilon of ee acrocs whe map should 
im @eneral represent vhe external variavion fiela om vhe 
Z=O plane due to currents Plowing in the lonosphere plus 
Tie wi Telas cue tO MOrmal IMbernal currPenvs wi tne wi enly — 
eoneucting mancvile. “~The vardation of Re across the map 
should represent in the main the magnetic field One whe 
Z=Q0 plane due to internal currents flowing in the Harth. 
Non=-linitormreres in these eddy currents: reitecy 1en— 


Unutormibies in the conducting Substrata. ~Quanvitavive 
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estimates, ofthe, depth and tondierivity of these non— 
uniform conductors can be made by use of the equations 


presenved inasecuions G4) and <1.5)s 


De eoeparavton or the magnetic fields 


peparaviol Of ra magnetic Tiel) anvoritesinverpal 
and external Darts Was first accomplusned cm ea eloval 
basis for the main geomagnetic field by Gauss (1839) 
and later for time-varying fields by Chapman and Whitehead 
(1923)3) batirisand Price (19390) Cand others. | ccparecions 
of Local magnetic frvelds on a plane have been) carried 
CUuG by clebert and Kertz (1957) and by Porathicucal 


ion O77 Oe: 


Global separation 


Global separation requires phe Tatvcine.of va maz— 
Hevic field known or assumed Over the Harta s suctace 
Gora auMber Of spherical Rarnoraes ands DNSwestameauzon 
Oconee Dut ONS vO each Narmonue Prom exvernal and 
inpernab sources. the mosy ceneral spherical’ harmonire 
representation of a scalar potential from which surface 
macnervie Taclds on a sphere cam be Gerivedwis wore tie 


form (Chapman and Bartels 1940) 
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n 
Ley 


wear 
Sa “tly itanhe 
eb 
c , VW « vial [ alae tai Tc sedaae 
re 7 
7 


a 
6.7) eastouet of beineeer 


} aa 
‘ 
pe 
4 
! < 
‘ 
ee 
- 
: pa uL 
‘ 
i y i iD) 
. 
>. A 
. Lhd ad! * i 7 
4 
tis ‘Vie i 1 > rig 
* 
an oe 
| j j j R i] ae 


(pi7aes:s avi e- ' / ear 


t a weer 


a Lie radi ussOoPT GCnee lero 


ia distance from the center of the Earth 
oa and aoe represeny Tractions of the harmonics oF 


cos mp and sin md between zero and one which are due 


Vee mSOUrCes OUvelde yhe Sphere. lf the radial gradient 
in tes 2 (tne vertical field measured at the surface) is 


Known 1t.can be expanded in spherical harmonies 


ay: ee 


r=a n=0 m=0 
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Bouiet 1 Onl >> seam be diftferengiated withwrecoect.tomim 


ia eee Uae Orme COMpDa inl cS OlmO te Ciidiom ice c(i Geain isi 
thee, Of equation, (1.50) will give the values of the 
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coefficients oe and Sake The coefficients A,” and 


a ue 


n ore found com, the Morthward, i) sand eastward, 


field Components, since these and not the potenvial 
COnotrtuge tChewobservyed data: 
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1 oV a m si tear 
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a m m eee 
= 2. x oe (cos6) (a, cos mp + 8 sin mo). 
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Y gives Une SOSPreLentva most weasily but & Comparison 


OF Ulcse derived irom A prevides@aleheck of Ghe analysis. 


Tals Mephodlis NOu apolacable tovamalysors Oo 
fecal ficids @madimly because it would require vod many 
Narmonies to-detérmine the potential accurately enough 
Overt =a Samrped “area . 

Price sand Wilkins (1963) made 4 three-dimensiona) 
PeUareavUgiwOL a surrace mMagneyic i Veladuby se of seuruace 
M1Getta le. we licseveny analysed Wwasstme Sait ielasoe 1932— 
boo S we eo horizontal Components Of the (Sq Bieta were 
tive mvOlaleo Weis a necessary CONnddt1on Ores One nulel 
ict mulaG cea bile iMbegitad sO Magicudce sel ce mamounGd 
any ClLOSsee Con voOuUlr One Une eiirlace Was Squaw yO Zere,, 
Tis Pelavron can De Tound Dy @applyime, sSvokes sine oren 
TO emer curl Tree Condi vom VY xen => wWeere fs the 
horizontal field. A grid was then drawn with 10° 
LAavicude “and iio EOWMeievUde WinwerseculOne sala yale smor 
Tinemtrebod. were found av lune grid pointe. | Once sa apo— 
tenuiewaWac C1Ven al anrbaurary valWve 20 -One poll a unc 
VOLpeDtad mateo Oylevspolmus Could besevalua: ed. 
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ie “soc OUMGrary Cone tary 

n = normal CO Une svirtace aus 

C= distance, from, pom Neto une surtace 
element cds 

» = ahy closed surface 

Vo o= Me a3 V; =the toval povential 


due VorVestine (1901). This Antegral was Used by Price 
and Wilkins to compute the difference between internal 
andrexvernal parvus Of Che povenvial atu any poinu fon 
URewsUrrace.. lh Che surrace "S as8a Cpmere equal omy leo) 
reduces vo 


Vo tweaZ 
ee ee ee 
i = saiweard normal componeny, Of force 


a= oaOnus Of, Lie sspmere 


peparavt on On va, plane 


The use of magnetic Variavuions Go determine upper 
Manes eoncuctaviIly SUrUChUre requires analysis ot local 
Pirekdse. Methods CF eanalysas Of these tields mave been 
and are being, developed. 

Siebert and Kertz (1957) attempted a separation 
ii eiwoo Camensi1ons cl the docal fields vector recorded au 


fifteen stations ecross the North German anomaly. By 
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findingscertain profiles alone which one componeny.of 
the horizontal field had no gradient, they were able 

to transform the three-dimensional problem into a two- 
dimensional one. (The separation inve aatvernalsand 
exvernal perts Was completed by means of an anvecral 
Operator. Sin many Cases This Condition is much woo 
Stricy vo allow a valid inverpretation Of the Separated 
mielde.) Whav its required tnen are anvecral formulae 
SpOroeriate for separation (on a plane in rie ieneral 
Case. 

Hartmann (1963) obtained suitable separation 
LoOriiikac Oye threavine Cie weparatsen Ol Gf welticeas, 2 
Drovlem an pOvenuial theory and Using Green's functions 
Of sone firs) ana second kind. Weaver (1963) developed 
two and three dimensional formulae by taking the two- 
dimensional ourier transtorm sof Laplace "Ss equatlom for 
the poventaal and applying ve Faltung theorem. 

Oldenburg Bion Poratnhn eo el. (1970)) have 
used the approach Of Price and Wilkins CIO Gs) irOmerrec.c 
a threé-dimensional separation analysis of local mag- 
Metic fields. Whe starting point for this separavion 
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For local fields over a few hundred kilometers 
the sphericity of the Earth can be neglected, imply- 
ing that the surface S$ is to be regarded as a hemisphere 


With an ainrInite radius, If (Xi); ve 0) Us an pont at 


iN? 
which the fields aré to bé separated then 


wee ean 2 5 ee FO nnd ame Om 
a A ape amelie fe an 57 
A etaken aS positive downwar - On the surface 3S, 2F=0 
and equation (1.58) reduces to 
- , + av 
25 Ol ae i Sn Oe pac Cro) 


DeirerenvrtatLon Or tie above equation will giverune 
Cope ee won or whew separaved magnetic field 2 vo 


give 


Ps 
N 


(X, iz Xi) . oT me ee 
P eee Re 1.60 
Coe se sl oe (1.60) 
Fate Oa OG 
(2, - Zio 7 25 Z ce ons 


When the equations (il JOO) <are combined with waco cbserved 
components X = Xe q Xs M4 i vo + Y; and 2 = Z ap Zs a 
complete separation of each of the magnetic components 


Ls obtained. 
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Geomagnetic depth sounding is concerned with 
Tour Gypes Ol Megnetic variations: 

I. Normal external field 

a. ANOMalous external field 

So. NOrmal internal rield 

bo Rou etionte internal field 
HU Ml-labicudes the inducing field is taken as the 
Sunol Che normal exvernal and normal Ginternal fields. 
pe nornalwinvernal TTela exnhibats Tae tadue von 
DrOpertties of Une “Tlar-eartm” models discussed, in 
Sseculons Wee and 123, while the anomalous anvernale t2eid 
fie co meOMe et Cra sEnnOmMorenel twee elie mahi etna Econ 
CUETO ac. Cbs cCuSced WM Seco On ah... 

it should be pointed uv that fields having 
SDatvial wavelenguns much Larger than the dimensions 
eT vine array cannot be Separaved by the formulae above. 
Thus at 16 the dimensions Of the array thavedicvate 
Witemeorre las are “anomalous and which are anseparabie 


Unormal' parts. 


iO iawn ve vons OL ceoMmarner ie denprl ys Oumd lie 


Ace WasepOllved "OUU ein Section 1. 9eaecomplevc 
Separavionsor the internal and external parvs Ol a 
transient magnetic field is not possible because the 


area in which observations are made is small compared 
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VO Tie Spatial extent of the primary inducing field: 
Therefore the anomalies must be interpreted to some 
extent on the basis of a preconceived mean conductivity 
GisvuributLon estimated {rom Ovner Sources. “The moce 
Commonly wsed source outside GDS as maenetotellurics. 
BUC “CONGuUCTIVity Cstimatves from magpnetovelilurics are 


greatly vstvortved am aréase of deep seated conductiviuy 


anomalies (Peeples 1969). This ts Gue to the distortion 
of the horizontal magnetic and electric components over 
areas Of exvensive laveral conductiviuy I2nhbomogenei ties; 
tie wOvane slavyered™approximavmon is ne, Vongverme valads. 

Poreavimey ah. (O70) haves stiowm tiation vine mares 
OL enVvesuI Gerd Onwds Laree inoue £oon estima veo canny tec 
Made Or wie COnGUCTIVILyY AaueGepel Prom separation 
resulUs. ouch GSvimaves: Gan be shown vO correlaverwicl 
Cceounysitcal models of the Upper manvle from over methods, 
mataiy seounermal studies and measuremenvs Om conduct ivi— 
TicssOn Uiltravascice Silvcaves 20 Nish pressures mand 
TeMmperaLvures.. 

im the real physical sivuavi0m Cvservarvons oF 
the meeievwe tmreld are made over suriace layers ihav 
Touma Malmmeconavicuing sneee (ol "real Compl e<iyy Beep e— 
Ssentingsa second limitation to GDS. The flow of shallow 


eddy currents may be irregular, becoming most severe in 
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tie socean Ghleci near coastlines and in areas of click 
scdimentary basins. The relevant calculations to 
G@ecermime whese effects are generally intractaple. —a 
peere, Nes been made by Praces (1940) in 2 study ion 2n— 
duction ain thin sheets. 

Figure (1.7) illustrates the method of Schmucker 
(1969) for a three layer model with a top layer of 
thickness GQ, a poorly conducting Antermediate layer 
Of Ghickness hand “a highly conducting mMali-soace,, which 
represents the Barth's suriace Vayers, htiehly resistive 
Zones Om tie Crus” and the wppernosG mantle eand ative 
higmily onducvuing mManvle respectively. ~The approxima— 
GLon O17 = (l + Typ, Cskin depth errect). Os =~ k 
(dependence on the dimensions of the source field), 
Kae eluate) paw weel Insbe Used aid ayers. 12 ces oc camiyic ys 
DeeteGo De tekenesnall compared vO tne wavelensun or ae 
Source Liela put Jarce with respect to the ekin-depth 
Valve pes thus the wave number k-of the amcideny i1eld 
€4am De ignored in the following relations. Ihe recurrence 


formule equation (l.2/)) gives 
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+ 67h, tanh Céid) 
86ih + tanh (6,4) 


Gao) 


with the resulting approximation to the attenuated CLEC Vee 


field at the bottom of the top Layer 
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fields. After Schmucker (1960), 
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The asymptotic behaviour of the curves in figure 1.7 
shows that p, should be at least three times the thick- 
Ness O fer che approximation im equation (ile 6m. tosnold: 
Mie wcueves DLOGred Tor Gitferane velucs. of d/m snow 
Chav veven) vhough pqe> d the racio d/m limivs, the uni 
LOCO ol fe ln) the bop ayer. 

A numerical example given by Schmucker (1964) 
nelps to silusurate the above argument. Given a4 kan. 


ocean (a, = 4x Le 


e.m.u.) above a highly conductive 
Menole sac cdepth lOO) im... then ds = 4) hme.) =n Oona 
ad/a =s725 48 sWhererore the atvenual lon of Ho aumogae ocean 
weliepe less then 42 10 pp > 12 kine eiving source faclds 
of about 38 minutes period as the lowest permissible 

Oe CLOd Or sri ce (Samer aod vero ld. 

Thespbase  Vvead ont the induced field to the induc- 
ing field can reveal roughly where most of the antvernal 
edaqy CUrrenus are with respec TO Uheim Ceol Inspec 
tionwot figures (1.7) and (1.3 ))shows thao af h/7d — i 
aeupstantiel pert of £G0 mand) HVO) are: alvenuated Tor 
values of pi/d < 1 for the out-of-phase (dashed) curves. 
This means that strong eddy currents are flowing in 
the highly conducting shallow top layer and vhat vhese 


€lrrents Will have a phase lead to H(0). As the valve 
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h/d and pi/d become larger the eddy.currents will flow 
On Che average ab greater depts: than dt. in the hMalT— 
space. [= 2s seen that these currents will produce 
frelas: shav are fia the main Am phase with HCO). | it 
Peale OOo sa bale. kom hen es bara yO Ord le t.d ie lidae ta 
roucnin between “shallow and deep biehly—-concuciting 
SuRUCDUreSs Simply by determining the phase Gditterences 


between the internal and external fields. 
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CHAPTER 2 


PREVIOUS GEOMAGNETIC DEPTH SOUNDING 
TN THE WESTERN UNITED STATES 


oclmucker (19604), Caner €t al. (1067 |) end Reitze! 
and Gough (1969) used Geomagnetic Depth Sounding (GDS) 
in investigations of the upper mantle Conductivity struc— 
SULey Ui we. COrd litera prom co Woof. Figure Cale wkocates 
the “stavions used by these investigators. 

Duringhone International Geophysacal Year iCiGy) 
1957/1958 the Coast and Geodetic Survey of the United 
pueaves Operaved 4 Net of seven temporary geomagne ric 
ODsehVaLOrles in Give wesvern United Staves. ochmucker 
(1964) reported that significant and consistent anoma- 
lies of the Z (vertical) variations were encountered, 
Scpectalinywauy Pelce Utam. "Norderarled study om tne 
Lyperia weCOnCUCULVaty anomalies Could ibe scarred ouy 
gue To tne @reauv Cisvance Detween vwne svations., “Wie 
average station spacing was about 400 km. 

During the years 1959/62 Schmucker carried out 
extensive observations av many Stations in the souvh— 
western United States. The Tine Or stations of most 
importance to the study of the Cordillera, shown in 
figure 2.2, consists of Six Stations across southern 


Arizona, New Mexico and west Texas. 
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ochmucker used the "base station" GDS method, 
with the permanent magnetic observatory at Tucson, 
me DZOnawais inis seis wor fixed Stayton, in is invec— 
tigetion~of whe Rio Grande anomaly (see section Wi7): 
He distinguished between the normal and anomalous 
batvs Ol vhe observed variation fields; definime the 
normal part as the sum of the smoothed external and 
incernal “parts Of Variations for the ‘casesot = ratner 
UninmerciiVe sternal source Tite Vd wand a horizonsaliy 
Seemed BCOnUUC TIN Ur yO MS tira PLO. wlaeum comma 1 
variations were taken to be uniform in mid-latitudes. 
DIncew nerve aS ami Near Comlrelepron veuWweem ves nomial 
and anomalous Parts, ochmucKker soughie vor derive ay each 
Seite OWvewCharacver Us UC sand: Rreduency-Gepen@eny sine ar 
Crane rormavilon Matrix, ©(P,@)),. Woe would express: au 
Cacimstaclon tviewanomalous ipart OF Une observed varia— 
hLons La cerms Of the mormal part Over the whole area. 
Hevderined the followine funcvLons Cover press. the cor — 
relation between normal and anomalous variations. 
PCPo oes Che time Varying disturbance vector of “the 
Cheer Gnvorietlom au le elac hon. 2 and a) and 


ie) represent its normal and anomalous parts: 
Ga) F(t) + Ce gee - Be Guy tr me Ce 
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It is assumed that a frequency analysis has been 
Deore cds On Vie O0SsScrVved Variations and that Dyn. Z 
are Nharmonue Lime Punctions of form exp Ciwt).) If 
ie assumed that a Lime Tunce pou, Z(G) as Jinearly 
oles LOS e CON mi Mme ei uneua ons Ce wa yl tbo oe 
being mormelmzed to zero mean value for a piven 
interval -%4 iv KOE RRs Ty: their relation will be 
independent. OT tine within Chis 2nterval and expres— 
Supike (oy 2 linear transter funccion Z(t) in che 


Peequency domain, assuming a eee ml Olen: Cie 


Ereansiorm Of Xt) can be expressed. by 
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Let S,, be the power spectrum of H(t) and S, the power 
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SDecUrUmNOL Gu) tien following wrelanilons, Bolan 


Sy = (Cy Cy I/T, 

¥ “ye 
Soy = Syz = (Cg-Cy a 
2 = Soy/Sy 


x 
where the product of Cy WicgimMebus Complex Con) ugaGe Ch 


x 
js ceal and positive and the product Cy with Cy cEyes 
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the cross-spectrum between Z(t) and H(t). The transfer 
Pune waren 254 is now expressed in terms of power and 


cross-spectra. Thus 
Z = - i ioe z(D)D. + z(H)H + NCE (COR 3) 


Where the (coetiicients in -equation (2.2) are complex 


and Of sone Torm 
Ze) Z f(D) 4 iZ_(D) 


The subscripts u and v are respectively in-phase and 
OUL-Ot—phase With the assoctaved dnducing componeny of 
Cie Occ wae Csan Or cine Wi dat ferences sila nae cuore 
if wewo-damensional field Gistributions equation (2.2) 


reduces to 
Z= Zo + ZX B)B + z(Z)2. 


with Zip =e (Db) cOosaut 40) sao (2) 
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Wheme B. 15 the projection of the horizonval disturbance 
O 
VECtLoOrEoOnue curectLon given by angle oa. For thewcase 


Orezero phase all quantities of equation (2-3) ere real, 


tana tee WH) Za) and for a 2-dimensional anomaly 
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oschmucker found that the Z(t) variation across 
the Rio Grande anomaly resembled D(t) and was nearly 
in phase with it. Therefore the anglé, a, can be set 


VOmZe Cored equavicmd (2.3) ives 
4 = Zo ot z(D)D. + 2(Z)2. 


Miewimelicavvon is thay one trendior tae monducitorm ae 
north-south. 

IieGecOnds “Or Geomacme vance substorm fields or 
Days Sschmucher Pound an increase in the 7/D ratio ease 
Ci COR and en anomalous 4 av COR Cileure 2.2 )— = the 
increased Z/D ratio at Sweetwater as well as the 
GCaducCedna/ DU rabtoty MicsoOnm are boun wore omsuidered 
"normal" in-the sense that both stations are assumed 


ie Les Wesiailigh Meee: Gugcehs Toi ImonmalaGiaure willy Fete ieeNecl T Oly 


Dire mie Menon COMO UCw TValt y) Sibiel Cr Ue one BUM eme Oo Clyic mean) 
Fevios sueecest. tndl either vne total conductivity for “ie 
surmace: layers 26 Consdvently Lower im souvnern) Wexas 
(tana soutiern rizonea or That the Uepthnwon tie oiehiy 
CGNGVic ie Waku et SUC IMO Der Mant let dks) eat em uineew 
Leta ee woen Uider Arizona. 

Schmucker interpreted the Rio Grande anomaly as 
Beevevemanic change Of bie internal parry Of tne 7 


Verte one Wiich es caused by 2 targe scale change in 
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Figure 2.5 


DenmuUuckKenes Second model Mor tthe 
Rio Grande anomaly (after Schmucker, 


1969). 
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the internal conductivity structure, since inspection 
of figure 2.2 shows that the higher Z variations are 
new restricted to a limited aréa; the 2 variabions 
do not return to normal as defined by the small Tuscon- 
bybe a Varvavion., Anomalous varcdavions as Geninea 10 
equation (2.2) occur mainly between LAC and COR. 
schmucker modeled the Rio Grande anomaly by 
employing the methods presented an section 125 toeevher 
with conformal mapping techniques applied to conductors 
OF ial ini ue. conductivity and of various shapes and) sizes. 
His firsc model. (Schmucker, L964): consisted ore noren— 
SOUGIeatoened SteD SUrUCLUre Wilin tne Cepia vo whe 
Mseetectesibsavitute conductor 160 Ln. West son UAC and 
320 kin. seast of LAC (figure 2-3) 1b can be seen) trom 
Pere ees etnata tne Sten model does nou explazae ihe 
sharp increase of the anomalous Z variations between 
BAC and COR, Which becomes greater for short period 
vVaruaplons wit @ More promounced megaul ve: anomaly 
between COR and LAC. The alvernavave model of tigure 
2.5 Gives much better fil vo the observations by placing 
the surface of the perfect conductor higher between LAC 
ana COR. Schmucker at first attributed the local anomaly 
AUC roe fighiy conductive., 100 Kn. wide iehaliow 


lever thatmucht imclude suriece layers, cxvendine norvh— 


ays trite 4 ctrsaty sits 4 
a 
: t 1d bs ‘dé pate) 36, a0 


j' baw iy gs pe with er 


aaa 


' 


Figure 2.4 


echmuecker's Ghird model Tor tne 
Rio Grande anomaly (Schmucker, 


1969). 
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south; between LAC and COR. A eloéser inspéction of 
fieUrevers SHOWS that it ts in fact he in-phase 
component of Zn which experiences the reversal 
between LAC and COR, which means that the reversal 
is) a2scoCiated with urrente: inva, bares, conductor in 
which eselt—-induction is more important than resdetive 
imowe biome. oechmucker (1969) has Gimce meyased fic model 
to concur somewhat with the above argument. Figure 2.4 
shows Nis New interpretation of the Rio Grande anomaly 
elones with CorrobpoOratinge heat flow results) along tie 
semen line as bis magnetic Stations. Schmucker Ss) anomaly 
Cemonstraves UNe=presence Of Surucgure In Upper mManuLle 
COnductiviuy av Une Dolndary Detween the Basan and Raneze 
anomunewarceay f laAinSs provinces. Tne 7enera 1 ormo1 1c 
‘conductive structure of figure 2.4 resembles the struc- 
ture (step plus superposed ridge) suggested by Porath, 
Oldenbtre and. Gough (1970) to accouny Por the Wasaven 
rons anomaly in Utah av the boundary between The Basin 
and Range Province and the Colorado Plateau (Chapter 4). 
@aner, Cannon @nd Mivingstone (1967) seu our an 
east-west profile of stations 750 km. long near latitude 
35 ° N from Horse Springs (HOR), New Mexico to Sayre 
(Say, Oklahoma (fieure 2.1). This protile was exvended 
aAnovner 250 km. eastward by Use of Une permanent obser— 


vatory at the University of Oklahoma, Norman. che 
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observations were made in an attempt to map the 
northward extension of the Rio Grande anomaly (after 
Schmucker 1964). 

Caner et al. (1967) used the "base station" method 
of GDS. Peak to peak amplitudes of recorded magnetic 
Verdetions an the range of periods 10 to 60 minutes 
were measured in three components to determine whether 
epsvatsen Was in a “low Z" cr “hngehe 7 zone. "Pine 
Classiiteavion Of the Stations in this manner was 


areed soy athe InGroduction of The ‘parameter 2. 


I = AzZ/[(AH)? + (AD)2]# 


WOLCM 25 The racho of the measured vertical To noriZzon— 
tale amour vudes.. The (ratios ay Various stavtons 
Gime red by FAaACrOrs Ob ONO OP UNres ande were. Gaus 
Casi ly Separated into cavegorics “Wow 2" and “nigh 1. 
Micwapproxamare Mocatwons Of Cher tow i lo eaves 1 
peansitvon Zones sien wane be readily wervermincd. 
Caner t»al. (1907) “used the wolar 
blOou method of Parkinson (1959, 1962) to define cor— 
relations between the amplitude of the vertical 
eempomeny and whe direction of “a horizontal change 
veCtore wit a hormal Svatzon the L-ratio 1s uncor— 


related with the azimuth of the horizontal change 
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Bast Front anomaly of Colorado 


(Gough and Reitzel, 1969). 
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VeCuUOr. 0 —arolan AD7 AN, and at an anomalous stauilon 
tie) t-raulo 1S a maximum Per Cercain values of o 

TnUs detiming 2 preterred direction. Ganer 4G 21201967) 
POuUnNC = nav Une Uransition sta vlen Sayre showed 
fevattos IMbermediate between those obvained ay ane 
western stations of Newkirk and Bushland and those 
aueuMe Measverm stations of Norman end Dallas. inc 
Stat OW oayre Was Cevord Of any detinive corre lavicK 
DeUNeene Une. Nora zontal change Vector ang une vert leat 
Change vecuor. A transition zone was Pocated) well vo 
the easy Ob the Rockies av che Vexas-Oklanoma border. 
Canter Piiverprered vliesemresU lips uac "."SWelle. WOU Memoacc 
Of the Rio’ Grande anomaly. Results of the 196/ array 
donot memo pore Gite wconcivsuon, | Later studies Porauia 
1969) SUS@ZeED thal (Periiam Sediments, Waien ara up) vo 
DPW teCcu witch In the area may Nave MeskeOmeae set teey 
On Une reer COndUcUIViILy Strucuure. 

Gough and Reitzel (1969) set out an east-west 
protile Norum of tna of Caner eb ally (1957) Vacea pour 
38° 30' N between 108°W and 104°wW, using variometers 
Of tubeir sown Gesien. Whis protile wasmoperavedsve 
search for a norvhward extension of the Rio Grande 
anomaly and was. placed with reference vo GY Stations 
eralecdvaliens( low-2) end Burlington (high-7) ri eurce 
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Figure 2,6 


Magnetovelluric measurement sites 
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Condiictivity vrotiles (Cancer 


Swift, 1960). 
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Summarized magnetotelluric 
conductivity-depth profiles 


NOTE: 

“INTERMEDIATE” merges with “CONTINENTAL” at 200 kms. 
‘CONDUCTIVE’ merges with “CONTINENTAL” at 600 kms. 
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(from Swift, 1967) 
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alone Chis profile. The H,D4% components of the dis— 
turbance are shown for CIM with the Z variation for 
pevenl stations. D is recorded in the opposite sense 
irom that used by Schmucker (1964). Z% as a function 
or cvime resembles D much more than Hs The transition 
eccurs at Salida which is about 100 km. west of the 
Pronuewnenge Of Une Soutnern Rockies. 

The anomaly at Salida was thought at the time 
UO pe, tne Monvaward- extension of the Rig Grande anomaly. 
However the results of the 1967 array (Reitzel et al. 
VOTO ROraun Ct.ral. 19705 this thesis) sugatest vias 
the ALO Grande anomaly is Continuous with the sWasatch 
Prony anomaly rather wlan that underlying .he =Soouthern 
Rockies: 

Swite (1967) ian the summers of 1965 and 1966 
Garrieds ou, large scale mazcnetovelluric measurements 
im the soutn—western United staves by combinange velluric 
data from seven stations with Tucson geomagnetic obser— 
VavOry doves  MVeuve ec some ves tie) LOCacuoms of etne 
Stawions.) Wie macnetovel uric dave were s0alyvecd stor, 
tensor apparent resistivities, principal directions 
and two dimensionality measures. The measured apparenv 
resistivities were interpreted in terms of inhomogeneous 
resistivity structure by employing theoretical values 


for two-dimensional models obtained with the aid of a 
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transmission line analogy due to T.R. Madden. The 
theoretical values took into account the known surface 
geology. 

The anisotropy of the measured apparent resis-— 
tivities was primarily due to inhomogeneities in surface 
conductivity and was corrected to give an areal pavvern 
that suggested an anomalously conductive upper mantle 
Deneavin southern Arizona and New Mexico. 

Mne conductivity profiles interpreted from the 
maegnevotvellunmic results were classified by swltio ac: 

"Continental" for Roswell 

iavermediave tor Pioenia. Gallus Demis 

MC onduc tive". Tor Gat ford 
andvcan be Seen in figure 2.7 alon= with the Ganewel— 
NeDonela prorile. “The Cantwell-MeDonala conducuivaty 
Vs Gepth prormle is the result of Cantwell (1960) “fitting 
Diss lone period magnevovelluric Wesily of coOnductivivy 
vs depth to those of McDonald (1957), obtained also from 
One Deri1od Magnetovelluric investigavrons,, and Lanimi 
ancarrice (1939). obtained= from the inc¢ertialt parvs. on 
semi-diurnal Sq-variations and smoothed storm-time Dee 
Varieuctons. 0b seems that the magnevovelluric evidences 
seneralily supports Schmucker's results in indicating 
increased conductivities in the upper mantle under the 
Basin and Range Province. | 

Swift has interpreted low frequency magnevoteliluric 
dataean terms of a Devrologiically valid upper mancvie con— 


duCrpivity Structure in a geologically anomalous region. 
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CHAPTER 3 
TECHNIQUES OF OBSERVATION AND INTERPRETATION 


Seo Le serray 


To The summer of 190/ sa program cb Suucy ofthe 
Conducvivity distribution in the upper mantle under the 
Western United States was begun in the hope of defining 
boundaries of regions of different mantle conductivity 
structure more accurately than had before been possible. 
To accomplish this, 42 magnetic variometers had been 
déestened “and constructed au the Sovuvhwest Cenver an 
Daitaswend avy the University ol Alberta im fdmonvon., 

Vie Aisi rimenins Were auraneed alone Tour Casu-Wesu ines 
avout ej 00 skin. = 10 Wenger hy Crossing the sou miernshockies 
and tne Colprado Plaveauy and extending weld anvo tne 
Great. Plains and Basin and Range provinces (figure 3.1). 
Though covering an area 700 by 1300 km. between 36° and 
4o° north latitude and 102° and 116° west longitude, the 
variometers were concentrated across the Southern Rockies 
ane wae. Waseuch Hault "zones The vevidence of verudcal 
displacements at the Wasatch Front and the abundant vol- 
Conve tor tne west Of 1G Mad eugtesved wae posstpit oy ol 
structure in the isotherms there. Schmucker (1964) had 


earlier reported anomalous or high Z@ (vertical) magnetic 


HAN OBS) a aeae alk 


The 1967 variometer array. 
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displacements at the IGY magnetic observatory at 
Price wsUGeal atl the Wasatch aul zone erea. 

jhe spacing between lines was about 150 lan. 
and that between stations averaged 120 km. Small 
COnCiuCtiMvi ty anomalies at depuis less than. 100) lam. 
Hieveve Massed py such an array, This riskwas 
accepted in order to search a maximum area for upper 


Menthe .COnOUct ivi ty Structure. 


See Varromevers and field proeceaures 


The variometers (Gough and Reitzel, 1967) are 
elassteal in type, Wilh magnets) (HD. 7 components) 
SUspemMOdcGeiOn TOrslom Wires. "Varlvatwvons 1m Shier maciet ic 
Pi eloe well oroduce Variable vorques On tnemmacneis wilen 
Cause Unem To move through small angles and the angular 
COs or OnsiOMm Wwhese, MAgneus she MeCOrced pMOvOcrapul easy 
inven lor Phorm at. ven second Gnuvervals.= Une mazmeuic 
Variavlons Will then appear as traces ol vecvangudar 
Specs ton 35 ml. tim. — TWwe Dase) lanes were Tecorded by 
meaneror ligne Spous reilecved, Brom nixed imirrors:. 

The most eeri1ous problem, Of 4 Classical type or 
Vaclometver 1s that of thermal response. This was Solved 
by encasing the variometer in a five foot aluminum tube 
and Dievying this am 2 vervical hole. 9 Calculations iave 
shown that a weekly variation of 20.6 Civ Gite os UG ecie 


O 
resules in a temperature change of only 0.10 av the 
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suspended magnets. A first order compensation for 

the consequent small changes of moment in the H and 

Z magnets was accomplished by use of auxiliary per 
manent magnets which were adjusted at some temperature 
vo cancel half the Earth's fieldrat the suspended 
magnet. 

The instruments were leveled with a Brunton 
compass to within O° Gn order to prevent contamina- 
GLOm ol TOne iconponent by another... Tae sencimiviny on 
Ciewitstrument ws DetTer Char jone camma andiwith @ 
Sample rave Ol vem seconds; relilaple records on 
MNe2oNevie FeVCiis With periods Sréeater mham one or wwe 
ibe se were Tobvained, 

Because tne energy "ol mapnetwe Varrvarvlons 
decreases gypeatly at periods below 20 seconds, the 
peooliem Ofvaliasing is Negligcible with sampling 
mcerve | OL ten seconds. “Theérmtiming was controlled 
by Sam Accutron timer made by Bullova, which wave contact 
Closures “every ven isecondsy every hour and ievery 24 
Noumea Small couls Near ihe mMagnere were sisedece 
derteceepne traces ao hourly and daily scontvacwmlosures- 

At periodic servicing visits (usually at ‘10 day 
intervals) to each variometer standardizing fields known 
to within one percent were applied by means of permanent 


magnets (two southern lines) or by means of Cols and 
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magnets (two northern lines). This procedure served 
to determine the sensitivities of the three suspended 
magnets and to give absolute time (WWV radio signals) 
which could later be transferred to an hour mark 
produced by the Accutron timer on the film record. 
From two such visits the rate of drift of the Accutron 
CoulG be devermined with the ascsumpuiom thay the abit. 
was “inear. “The Accutrons were found to have reasonably 
COnstanv raves Of order 10 seconds per day or less end 
Lime Was @eneraliy relivable within one minute. "Wich 
S00G iweld operation Detter than UwoO percent preci ciom 
of calibration Can be Secured Dy these metiods. 2lticough 
Diem ole Cision auuouned av sOmMew aval Lous Wace sles ce: iam 
Go y 

Anmetpdependent check et the relavave: time pet ween 
Phe .nsvrumenus Of bhe array and also a check vom the 
assumpuron Unease the rate of driit of the Accurcon was 
approximately linear was obtained by correlavion of 
geomagnetic micropulsations which should occur simul- 
Sanecuciy (Within a Minuve) over the entice array 
(Rostoker, 1968). These micropulsations, known as 
Pi2's, have periods 40 to 150 seconds and are clearly 
visable ab the beginning of ach magnevic event yay all 


Seat rone on che NorwZOnual ULraces. Thus the maxamum 
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error in timing was mo more than one minute and sin 


most Gases Considerably less. 


5.5 shlactielLc Variation. recoras 


This thesis includes analyses of two magnetic 
variation events: 1) a magnetic substorm occurring 
from 0500 to 03830 UT on September 1, 1967 and 2) a 
moderave storm occurring trom 2300 to 1230 UT ‘on 
SPepremper 20-21, 190/. Whe substorm was selected 
BecaviccnoOr pune large amp ludes Of Marner Cuvaris pron 
ano, because Ll Was recorded with nearly GSZ.etticiency. 
meso lO Wasuselecved because Of shUs Menevheand vworied 
ireeQueitiey CONventl, DOUN Of WhiCciare Cestrabie TO. Demme 
pasmiserOlr sLOne period Varivatvaon navelds Sin Scud yams 
Secuc hires = line ricrency OL recording, Tom ume: svorm 
yeried Lromveast: towwest. | “lhe western halt toh the 
aeray . shat covering vhe Wasavch fault Zone wrecorded 
Wien Bop, -Or preatver erliciency and Une veas tenn hatin: 
that covering the Southern Rockies, with 63% or greater 
efficiency. Greater efficiency was achieved in che 
horizontal components, mainly because the Z component 
was metastable at several stations. Results have shown 
Gham the Horizontal components yield as much antormatvion 
abou. the structure of the dnotialies as the wervical 


component, particularly after separation. 
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Endarged Xerox prints of the film records 
containing the magnetic event were comercial ly 
digitized at one minute intervals. Before digitizing 
the Xerox traces were smoothed by underlining them 
tnepencilsse as to avoid aliasing trom micropulsataons 
Ol Perlocs! less than abouvetwo: minuves. eine amar a0 
value of each component was subtracted from each 
Succeed Me value. Physically thie means vhav the 
freld is thought to be quiet before the variation 
ana the initwal value can be sev to zero. The com-— 
ponenws Of the variation fields were normalized by means 
Giviacvors calctitlared trom Une field Sstandardiza vem 


data. 


3:4 spectral analysis 


Geomagnetic depth sounding as concerned waivun 
Variations 2 Conauctivity eb depuh. Consequently 
the skin-depth, p = Wiuaank. becomes an important 
paramever to estimave. The penetration of vine external 
field anvo tne Barun, the ancue rion process iat con- 
ducvive structure and the €scape of the field sot the 
impernal cureenve depend On vine Trequency,,.-W.s Tie 
relation between the external and internal fields, 


whose vector sum forms the recorded variation field, 
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is likewise a function of frequency. More information 
can thus be extracted from amplitude and phase rela- 
tions at definite frequencies than can be obtained 

POM Ne Orieine | (daca im ohel Game domain. Hor these 
reasons it 1s desirable to represent the total magnetic 
variation field as a Fourier series in the frequency 
dome im-ana vO use the Sine and cosine cocriicicnts of 
Ulemwscrics OF Une amplitude and phase, to Gisolay 
Spatlal Variatrions in the mapnetic Variapion t1eld 
ecroOscs  oHe wrray. 


Maks obo aeNa" 


dt Ge a8) 


Peau Neeeene ral. Vormulavionwor FOUrmersiturans hOrmm Ol ca 
time series/X(t). The Fourier analysis used throughout 
tis ness. was: bDhavasiaven by ihe REARM program gin, one 
IBM SSR package employing the Cooley-Tukey algorithm. 
Given 2N real numbers Xo. Xs» X5 ates Xoy_y ths Ssube 


routine computes Fourier Sine and cosine cochaicicnus 
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Each magnetic variation, which can be treated 
as @ time series, had a linear trend removed before 
fourteen analysis. ~The first and Vast ten minutes of 
Sach vimel Series were smoothed by multiplying them by 
COeiPiclenvs Of 4 Sine function Gvaltateda from 0 to 
We Vint s smoothing is necessary vo prevent pie: antro- 
Guction of Tictitious high frequency components that 
would be caused by sudden changes at the beginning and 
SNdgon “Aeiiniver time series. 

Ther Past Fourier alcorithm demands that. che 
numbers Ol points in the bime Series’ be am antverral 
DOWer Or twoe 9Thus at) wes necessary Go add Zeros. co 
Tice DUT Ommvne ame serme>s —vepResonG Mie rMc ws lommo ln 
andgcunews torts. IWurunier Zeros, Deyond Une siexoy invezrad 
Dower of Two above the number of dava points, were added 
tOrnerecase the number Cf Spectral esvimaces. Brhysi— 
ce llyotnis te reasonable Since the subsvorm = and Storm 
can be regarded as transients in an otherwise steady 
Macnhepiec Tield, so thau whe sveady magnevie Tield.is 
represented by a zero amplitude time series. 

Toerind tine effece of “adding zeros; en analyses 
was eae ee Dv Mir. Daw. Oldenbure and the writer in 
whieh the length of the work vector M (the number of 
tealsdaca points N, plus the number vol zeros added) 


took values 256,512,1024 and 2040. N was 211 for the 
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substorm. All frequencies found in the Fourier analysis 
when N=256 were found for each succeeding vyaluc or M. 
It was found that the sine and cosine coefficients of 
the Fourier series were attenuated by different (for 
each ditrerent value of M) but constant factors. “For 
the coefficients to approximate field values in gammas, 
each ycoeiiicieny musy De multsa plied by a sfactor N/M. 

Ti vais eis done the sine “and cosine coecttictenne wil! 
be vuhe “Same at. 2a particular frequency for ale values of 
Mae Whive addivional esuimaves secured by aAncreascang iM 
(ust lewwess. relaable. Une scectra or valves Onl = ticom 
256 to the maximum value 2048 remained smooth. For the 
present analysis M=2048 was used for the substorm and 
M=4096 was used for the storm which is represented by a 
time series N=81l. 

Mebane doLy Vora eneck sone uMe Vai dik yy eon emcee 
of the Cooley-Tukey Fast Fourier algorithm with numerous 
zeros added after the real data, data from the substorm 
Me Renee i, LOOT recorded ab Several esvecionsmierc 
Fourier analysed using bovum tuhe Cocley—Tukey algorithm 
and =a slow Fourier program which did not place special 
restrictions on the number of Gata poinvs analysed; 212 
points were used for the substorm without zeros. The 


two methods yielded identical spectra in both the 


23 


; | : 
: | 7 7 
fe 7 7 


_ 


; 2 —T ; an mG j : 7 ad fs ie, iD <a AY 
- 
| a | 2a or orn 


a, § fumed tow rl 
* 7 eo 
yrne @@lves. adhe aT enc 


oy Gorbay sewseT bap Ss 


mf alas sas 


a . ; atusiotttsoo eae 
, : 7 _ - 
j 7 ol wel ORS > 7 


Pa oure 3 2 


Fourier spectral amplitude and phase 
hor Very. cal Componeny Of subs corm 
of September i, M967 {8 tor starvions 


across the Wasatch Hrone on tline i 
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POUrLer  soeci ral anlolitudesandmpaase 
ior Veruca l VCOnmponen) © asUvs uormm 
of September 1, 1967, for stations 
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amplitude and phase. This check was carried out by 
Dr. H. Porath at the Southwest Center for Advanced 
eludes. 

At some stations where the amplitudes were Very 
laree  Meiniy on lanes 1 and 2 for the De component. 
fictitious values were interpolated in the gaps seen 
inthe Calcomp plows in figure 422 to wermic some 
evvempt at speceral analysis. The results for these 
stations at which the gaps are numerous should be treated 
Wit reseryc. 

Bigures, 322, 3.c8and 324 Sire examples of sine 
FOlrler vranss Oris Ole Une Veruleal Componentes ol sume 


SuoSstorm and che siorm form some stattons Over the Wasaten 


Prong. | The substorm Nas. a. smOoouh amplivude specurumn and 


VErVeMearked [Osi Trerences Aare ObDSserved) between, svations 
PORN pe eLOOs SreOMecOUMinvves BuO 2anOurs. s NOs trope Cure 
Contain very Ligtle energy au periods Tess  unan 20 
Minapes. souls tanctial phases dit Perences are observed 
between stations, especially where anomalous induction 
Occurs. Bub Gor periods exceedin= 90 minuves, Unese 
phase differences may be partially due to variations 
itm tae one period trend removed berore Hourvernsanalycis. 
A tine specorum Giigure 3.1) ls obtalimed for ene 


storm with peaks at periods of 155,120,85,65 minutes and 
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Figure 3.4 


frourLer spectral snp ligude sand onase 
Lor Vertical COomoCnent ~Ob isvOrimmro 
September 20-21, 1967, at stations 
Of Hbine Ss on eather side of the 
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several peaks below 60 minutes. The spectrum again 
Contains very Tittle eneresy for pericds under 20 
minutes. A differential shift between the energy 
Miexiniawor the Dorlzontal and yertacel components is 


observed in the storm spectra for periods longer than 


pO MMinubes tigre 3.5), Thus’ is probably telared to 
Cet eles anduced by Anternal currents, | Mie wim terias 
Pleboe becomes, relavively smaller Gor loneer period. 
Quey To one ower eriiciency of the Gdncductwon process 
au. tne longer periods. Above a conductive halt—space 
Divine tive COndUCULyLLy Varles Witt depch alone. tie 
NOrs Zone alate lds (Of ihe eriale Clm icetl (CmalLiicCieas cai imc 
external horizontal fields, wWhereass for the verpical 
component) ne tnvernal ile ldeis (Opposed vest Ne se xtemial 
too timawdacne Lresubtansg cCompOonchieans thetdliterence 
(Cselcvilcn le. = “Une neu entecy is a shvrp on the peaks 
Ch te Morizonual cCoMmoonents woward tle storver periods 
anc conversely the vertical peaks shire toward une 
onger periods. 

Regular varlatvions an amplivude and phase of 
thewweepectra av Neighbouring Stavions Ssugees vy] chau 
EPiesemenavial warlavlons may rellect large scate 


conductivity changes across the array. Results of 
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mapping these spatial variations for several periods 


are shown and interpreted in chapter 4. 


53> ecOeration of the variation fields 


Povevaluave the surtace Aenean Gn yequat: Lon 
(55>) quoted, in chapter 1) the magnetic components 
Must De pl Oeved “Ona rectangular grado) The stations 
wece -Uransierred from a sphere ¢O a2 plane by plocr ime 
them on an equal area map. The Z-direction is taken 
normal eto any porn on whe plane and the % and. ¥ 
Garectlons are vnose ot the lonearuce snd wacictude 
auevtave point. “As ube meridians are mov parallel. 
aS ceansvortiaulon ©O a Cectangalar srid ms reduaredr 
Mane wseparaviton OF the- total Mapnetie Tvela Variataon 
LivOmtrS im erie and external DerGs, VS CMs Carried 
Out on the new plane. A north-south reference Line was 
drain Birouch he center Of Ghe avray, and ay eacn 
Svat on une enoele nN bebieent a dine para lle lto sence 
fererence line, and the meridian indicating, true geo- 
erapnic NOLr a. Wes sLOuUnC 92 iss Gn Variant Sor tai 
Granctormation. Let XP ana YP denove northward and 
eastward field yalues on @ rectangzuler grid parallel 
to the reference line. The transformations vo convert 
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x sin 7 = Y¥ cos 7 


ihe transformation (3.3) is applied with aopre— 
Oriave fH vO. Une Tame series at each station. “The 
Uranstormed vame series are then Hourier analysed and 
mapped. Por tse remainder of This sect lon co avord 
COntusron x. will denote XP and Y, VP. 

TO draw Concour maps of the sine and cosine 
transtorms of the X,Y,Z2 components, the values between 
Sta cLons muUsL Ibe interpolated. “As @ Lirse appeoximarvion 
Une wivekds: are contoured Using Minear inverpolaraon 
peuween Svavlons. A second Constrainuy cam be applied 
COe Une UN berDOlaGlOn, DrOCeSn a 1) 2S) requingeds vial etite 
fields be d¢grivable from a scalar potential. 

Diapie wield Ts andecd Werivanle rom a scalar 
potential function, then the magnetic field H = H(X,Y,Z) 


measureo. on the Barth's surface must be curl—free ace. 


Vx BS 0 
Therefore 
COL oY nO _ BWA 1 ao J OX Bi 4 


Bech term in equation (3.4) must be equel VG Zerc. 
Because observations were made on a plane,only the 
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ax oy (83.5) 


Wal provide a ywsoable constraint. Therefore aruer, che 
maps of the XP and YP components have been contoured 
Dye tineear 1nverpolation, Tnescontours. cCancbe ma juseced 


TO me Cue tierCUrl= Tree ICondirion, an eCwabwolu so. 
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where 
Le = distance in the x direction between Y(1) and Y(2), 
ly = distance in the y direction between X(1) and X(2), 
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Practice has shown that when superposing the 
Manag -Y maps the condition (€3.5) was satistied best 
when the gradient of Y in the x direcctiom was large 
(Ofcourse this means that the @radient of < dmethe vy 
direction must also be large). If the gradients were 
emai tie eur) —tree *condiuon gaye Mitile i one way 
Cf ys mprovemenl to the incarly ianterpoleved sconvours. 
Miverpolariton On a2 Convoured map using fae Clg ieee 
Condition was first Ssuegesved by Reivzel (Gougnm and 
Reitzel, 1969). 

pince the CONUOUrING Was Carried Out On am x—y 
plane, $= and a are MO Pacer ance NG sa Miata te 
eonmevreaint vo equation (3.5) ean be employed aim che 
COnmvOuUrine Of 4. 4 maps prepared Tor separavion were 
Conromred py Vinear ViverpO lated On ane modi lode py, 
avording Sharp changes in gradienc. 

The =magnevle fvelda Componeny Maps were exvrapo— 
fated beyond tune Limits or the array to give Detter 
Sstimeores or wiervalics Of Vie sSUurTace Gnvesrals= aL 
Pomme Wetiim the array. Whe curl [ree Consvreiay 
Wae seo on the extrapolaved horizonval fields. (such 
an extrapolated field will clearly convainm errors and 


Psecervaimiy mov Unique, but should give a Devcer 


A0prosimacvuon vo the real iteld than, tor insvance, 
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the assumption that the field was zero or constant 
outside the array. The extrapolation was carried 
vo 300k.“ outside the array. The farther one extra— 
polates the more unreliable the field DUT. The=errecr 
of errors: in distant elements. of area falleeas aly ae 
(His the distance from tne’ point of separation to 
UieeeUriace -clemeny, equavion  Cle5>) eo ihe 4 tela 
Was extrapolated by smoothly extending the contours 
CutWard trom the vorrmeinal “array. 

ine surtace integrals may be determined 17 che 
average values Of “%,Y,2,0%/0%, eY¥/oy, 302/0% and 8 Way ere 
known for each square on a SYMMeEUric Grid, Centered at the 
POlNU Tor Seharation and used as a “summavion’ window. 
Oldenburg (1969) evaluated the effects of square 
Windows Of sides (00, 900 and 1100 kms; and, found 
Caste wne=resutte even wiear! thes ecdgees, Cle therarvay 
var ved. only by) 1 vor) 2 gammas. “Oldenburg (1969)) alse 
Wreovel 2 program to compute the suriace invegrais by 
Moving ea -square rid abouv the maps.) Tae valuesor wie 
Gifrerence between the internal and external parts of 
each magnetic component was then found for the center 
of thie “window". Near the edge of the array the 
separavion as Gependenl oniahe lessyreliable exvra- 
polated field and some values inside the array are 


Srolmdedo. in this way Cons lerab Ue ernore sat bec mlea. 


the edges of vhe array. 
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When separating in the frequency domain the 
eosine end! sihe transform maps musi bé plotted, 
contoured and separated independently since the 
Sublace anveprals can be applied only to fields in 
Gonsvant phase across the map. After separation 


yney May- pe combined according vo the relations 


: F - Sua ie 
= 7 2. ee ead 
we eos oe + sin Pe) ae ea e a ne F 
(3.6) 
. . L Sala Be 
= + i 2s SSS SS 
es (cos a sin"F,)* ; tan Deters nas F, 


where F can represent any magnetic component (X,Y ,4) 
and wie subscrapus © and A Teter to the "exvernal and 
tiverialeparrcoerespeoulvely.. 1G snould pe sien done 
Uhatewhe numerical magnitudes or Lhe cosine wand sine 
transtorms have physical meaning only an relation to 
the sien and magnivude of the spectral Compomeny 
Concerned at tune stare Of Phe wWime angervaly Jal 
Eranerormea. A change in une srarving.time will 
invroduce changes in the relative magnitudes and 

Signs Of the "COsine and Sine vransioOrms, Vous devermiina ng 
whether a Piven term in the time series appears in he 


Cosine sna/or Sine wransilorm or vin “Dotn. 
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CHAPTER 4 


ANOMALIES IN THE VARIATION FIELDS AND CONDUCTIVE STRUCTURES 


4.1 Interpretation procedures 


Three components of the magnetic vyariavion fields, 
plotted by the Calcomp as standard type variograms, were 
first examined for amplitude and phase changes along each 
casu-Weey line; especially in 424. Next, maps of three tyoes 
Were Dlovted Tor Cditieren’ purposes: 

(a) The amplitude and phase of the Fourier transforms at 
selecved perlods were pilovved and contoured.) These maos 

Of the unseparated variation fields have been found very 
eonvenient for locarvon and Semi-quantitgavlve aniverprera— 
UOMO en verial ICONdUe tye Tour uicuUre s:. 

(b) The sine and cosine COctMIC Tents (OR Une  tOuUrter 

ULans nore wou weelOCvuCad ™ requenCc bes Were Mapped. and sConvoured 
BOG sera cavron On Whe Mvelds 2nlo exvernaweaad TiGeraa 
ete see er Od domed nm. 


Ue) inetenvaneous tleld components at selected vimes were 


mapped and convoured by DW. Oldenbure (1969) for separa- 
Pcie iar ou limes COne lia, 

The procedure developed and used by Oldenburg (1969) 
Nas been used to separate the Magnetic Pields i1nid internal 
anidwexternal paros; (bhe invegral expressions used, are 


quoted in Chapter i) Coquations 100) 6 Mor ssecara toma: 
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the time domain the procedure is very sensitive to 
GLscrepencies Detween Stations im tame snd the times 
must be chosen when all three components have reasonable 
empurvudes. Ac least four ditterent Gimes have co be 
Separaved CO BPive any picture of he internal ane exper— 
nal variations of a substorm; Oldenburg (1969) separated 
Che Subsborm of September 1, 1967 at the times 0030, 0645, 
OPO} (a> Ul, The substorm was Tour and ome—nait waours 
long and it would be desirable to separate at more times. 
For a storm which is approximately three times as long 
eng sor Varwaole frequency Convenu a pronuprraye mumoer 
On InNsivanvaneous fields would haye vo be separeved.. As 
each SeveravioOn Ol Ansvanvaneous twelds involves. Gace hand 
CONLOUrane Or Six Maps, Separarton au Many Tames 2s mo. 
precu@eanle by Our present Meuhods. (separation in une 
Gime domain gives a very Pood picture Of the Gime vyariacions 
imetne Separaved parts across tne array. Ine exvernal 
fielavcomponents provide esvimaves of the scale Vengun > Nn; 
Of the external field. (Results also illuminave the mor— 
phology of the auroral magnetic event which from our 
Viewoolnts 16 the source function of the input Signal vo 
the array (Oldenburg 1969). 

jhe above separation procedure 1s also valid ain 
une ~“Erequency OL period domain. The advantages of a 


frequency domain separation are: 
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(1) Time discrepancies between stations are less critical 
chan ii tae Gime domain. 
(2) Phase and amplitude representation of the external and 
imbernal fields are Obtained at the chosen periods. Thais 
is important because magnetic mane and. Gill tis ion 
processes are frequency dependent. 
(208 The Fourier transtorm maps of type (>) show tae 
Conductive structures more strongly than the time domain 
maps (¢), because in effect the energy of the whole event 
avscach selected Trequcncy 1s garnered An che Fourter 
Dransnormavion., “—POr thas weason the discuss von on mode is 
Cie COMdUCTIVe SUrUeCtInes makes Se (Ol Scparaved elds na 
tie (erlod domain. 

The Use Leone VarLous Dypes Ouranoma yams, 
Sepabavlon O01 Maelds and modelling oF conductive setruc— 
tures have been discussed by Retuzel, Gough, Porat and 


Anderson (1970) and my Poravh, Oldenbure and Gough C1970). 


4,2 Variograms 

Variloerams or a2 SubsvOrm (o% Toeprember = 15 1967 and 
a storm on September 20-21, 1967 are shown in figures 
Wea, 42a and 4.25. The east-west lines of variometers 
are numbered 1 to 4 from north to south (figure 3.1). 
The variation anomaly associated with the eastern front 


of the Cordillera, formed in the area of the array by the 
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Gaede go Variograms for substorm of September 
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igure 4.248 Veariograms Or SvOMmM On Ipepvemoer 
20-21, 1967. Two northern lines 


Of eaieice vee 


ba 
1 iy! a0 es okie ae 
“co A iglesia ¢' 

| | ; iV 1 4 _ 


ae 


tL SNia 


¢ ANI] 


L961 "Oe Y3IGWSLdIS 


de 


aaMats 


yael OS F 


t 3v1J 


SPREREEE 


amalal 


Figure 4.2b Variograms for storm of September 
20-21, 1967. Two southern lines 
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southern Rockies, can be seen in the vertical component 
of the substorm in the increase of Z between CAS and 
CRW (line 1), between LYO and STE (line 2) and between 
SAG and WLS (line 3). Only a small increase is observed 
fom the substormeat theveastern mosh station, TUC mim 
line 4. The presence of the East Front anomaly is more 
obvious for the storm. This difference in the response 
OPetke Vervulcakuhieldsror che sborm as mosteuce lydue se 
hive sent ereny aZzimuun Of: are external Thelds: the dominant 
azimuths are NW-SE for the substorm and NE-SW for the 
seorm. "A further bay event on September 26, 1907 also 
shows the presence of the East Front anomaly in northern 
New Mexuco by a substantial ancrease ian 4 bpevween) UVG tend 
TUG (Porath. private communication): 

tie wWVarvroOsrans ws trolly SnGiCavew ble Bo xistwenc emo 
another variation anomaly with roughly a north-south trend 
6Verrcne Wasatch Pauly Zone in Uren.” Dns wall we rmalived 
thesvasetch Front. anomaly. [he ancrease in the 2 component 
can be seen between the stations MUR and MON (line 1), 
Piipoanc MRA (line 2), SUN end EMe (line 3) and VEY and 
MCo wcline 4) for the substorm. “The Z-component Cf the 
aeloml Waa not recorded at HMES Duy increases im 4 can pe 
Heen petween the Ovhner pairs of stavions. he similarity 
in the traces of Z(t) and D(G)y especially dn migure 4 ys, 


indicates the presence of anomalous internal Z at these 
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Stations. Upward (negative) Z variations are seen to 
be approximately in phase with westerly (negative D) 
horizcntal variavions. This and the fact Ghat =tacions 
west of the Wasatch Front (PIL, EUR, HIK) have very low 
Z amplitudes suggest that the Basin and Range province 
Usa Lee von of nigh mantle conductivity. The Calcomp 
plot of the storm shows reversals just west of the Wasatch 
fronvu, | [ne -piase reversals “dni helZ scomponeny tor per ods 
Of about 30 minutes av BAK and PIL, Sugeest that “a ridge 
ia che conductivity structure lies alone the edee of the 
Bapiieeanc hange province at the Wasacchm Mrout. 

BRU at the west end of line 1 shows Z(t) resembling 
@ horizontal Component with contributions from poth H(z) 
and D(t). The variograms suggest upward Z following 
SOutnward Hi. This may andicate the presence of a highly 
COndveving body south Of this stabion, which, could pe 
inverpreved 42s the morvhern edge “of vne Bacin and Range 


prov Lace. 


W3 Maps of unseparaved spectral Components 


Figures 4.8 co 4.6 represent Maps of Thesspeetral 
components of the substorm. Figures 4.7 and 4.8 represent 
Sitiular maps Of Une =svOrm, 

The vertical field amplitude maps of figure 4.3 


show two prominent north-south striking maxima, one along 
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Figure 4.3 Fourier spectral amplitudes of the 


Veruca ll (eomponent subs borm fiend 
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Figure 4.4 Fourier spectral phases of the 
Venu Leal COmponeny slvacolm 
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horazontal components, of Ge 
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Figure 4.6 


Fourier spectral phases of horizontal 


€OmpOMents Ol the Subsvorm trelas. 
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the east side of the Southern Rockies and the other east 
of the Wasatch Front in Utah (see figure 4.9). These 
anomalies show more strongly at the longer periods partly 
because of the greater penetration of the long-period 
fields, but partiy also because the magnetic events have 
Moce enerey av the lorger perieds. “hab the second ~nrect 
De Mere amportant in the present case Ss sugpestea by the 
Virtual disappearance of the anomalies from @ man of 7 
amplitude at period 120 minutes for the substorm (figure 
ie 0} The storm has high enerey am the Z spectrum near 
pertod J50 minutes (fieure 3.4 ) and whe anomalies show 
well ine the 7 amplitude map ab this period (ficure a7). 
Les Iimporcany WlO select Uertods am Which tie selerey ais 
Drei. som whe Sepacavivon Of wields and model ycalculara once 
Garermust be taken) in tne Use of the Sscorm maps 
to remember an inequality in the amount of information 
available Tor contouring the eCasvern and western halves 
Ofethe array. About ninety percent eiiiclency an recording 
Was attained for the west and only sixty percent in the 
Past. siliso at several of the stations in the N-f secror 
GP the array the 4 and D vraces were driven olr scale by 
the large storm fields. The gaps in the discontinuous 
Ceaces OL 11gure lee were Billed with Tictiptous, anGger— 
polated data to give some attempt at spectral estimates. 


The north-eastern parts of the maps of Z and Y in figures 
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Figure 4.9 
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pimplifiedemap or two local internal 
Currents and welavionships, peuween 
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Figure 4.10 PoUrier spectral amp lriude ore vertacar 


Component (Of SuDStorm avy L270 minutes. 
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4.7 and 4.8 are therefore to be treated with caution. 
Greater importance its to be attributed to the results 
ebGa med irom une western Half. «Since thie 7s an ares 
700 km by 650 km the results remain important with 
Pecpect, Co deve mining Widerlyine Decent SULUCGUREe.. 
Tme Mapse,Ot figures 4.3 and 427 tor 4% acy T=45 minutes 
for the substorm and storm are strongly alike, suggesting 
that the contouring was reasonably accurate. 

ive “Casctward Horizontal Component YY shows sprone 
evidence of both the Hast Fronb and. the Wasatch Bront 
Cupeent = svysuvcis dma thie amplitude maps.20r cae suostorm 
(eieure 4.5). She maximum Of bY welwesteon the been aaa 
DO@mIMwCas Cone Sil TCaClu Ca sem Wei Ill Sata hile On ule mine. (lao iii 
Maximum, Indicating a linear current. Al Ene Waser chy iiront 
bole suo DelleVved TO represen: a local upwelling Jae ene 
Conduct uve Manvule i a widee Superimposed on ey svep watch 
leaves une Conductive medium Higher under the Basin and 
Henge provitice than under the Colorado Plaveau. | Schmiciker 


GhOGO Nea chned similar cone lisions walt respecu TO Uile 


Conuuctaviuy at the edge Of the Basin and henge soroy ince 
in southern New Mexico (figure 2.4 ). Later sections will 
Some that similar GCepubs 40.5 uae) Upwelling were falsou ound, 
suggesting that the Wasatch Front structure is a northward 
expenciom of the Rio Grande structure. Tne Y semplicude 


maps at the East Front suggest a SUP OL Clare nb aUnG els 


pe bry 
me | | 
me | 
a 
* a - ' : 
| Dy y OAé a aro } a 
in i Atnyatee Se I qa9" 
ign? $i ct Sd 
© _ 
= 4 ! wr sh, 
by ay Seat 
lp J 0h unputiedimc et? 
= = By 7 a 
py ay 
‘Lar 
iG 


ie af 
A { i - btie > ( " Tepe 5 
op) Meiged iso ware av enaw cate dad 


- 

7 fu) mn 4) a ‘hekaT- Je tyr) 
io? ee ie ree nw 
; , we! —_" 
nay 


12 


124 


the Southern Rockies between the East Front and the 
Colorado Plateau. A comparison of the Wasatch Front 
and Hast Front Y anomalies at T=30 and 89 minutes shows 
Chats Che width tor the esti p woulun have vo be wader junder 
tne Hast Front than under the Wasatem Front (fieure: 4.5). 
The amplitude maps of the X component of the 
substorm (figure 4.5) are relatively featureless, as 
would be expected where the major @onductive structures 
StUrine NOrthn—souvh. Wi ismall anomaly umm Oe associamed 


either with the Uinta uplift or with the deep sedimentary 


Uinivea Basin 20,000" ft.) am ortherm Uva aie tscenmiis o 
sine transform map at period 50 minutes (Reitzel et al. 
1970). As Schmucker (1964) \ocints out, deep basins: of 
Conducving Sedimenvs Can 'cause anomalous Varitavlonsmat 
periods up to one hour or more. Superficial anomalies 
usually show a Conspicuous phase lead of the induced to 
inducing field Rcenecten: 1964) and early arrival times 
ors for the Ssubstorm in this region may be an vexpress pon 
Of the Uinta Basin anomaly - 

Tie X% eomponens maps Porstney storm saow 2 reversal 
in north-south gradient of the amplitude between T=45 and 
150 minutes. Whereas X diminishes southward at period 
Woeninutes 16 inereases southward over the entire array 
2G t= (50eminuces (iigure i je SOpher mapper auotmes town > 


4ndicate that the reversal occurs close to T=90 minutes. 
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By way of analogy with the substorm it would appear that 
the external field should decrease southward. It must 

be pointed out clearly though that the analogy may not 

be valid in the absence of 4 model of thie Ss Vorm Veurrent 
system. Discussion on this feature of the storm transform 
Mapes Wii! be continued in the mext: section, seieamedealing 
With the Sséparated storm fields. 

Une Diasec Maps NO only conraimn inmormavaom wan 
regard ve anomalous Conductivity Structure, DUT also 
enable tne sieniticance Of the amplicude manos to be 
eeoes ocd. wena pliases Ole and Y vary smoounly acrocs 
Sec milapws A svrone Z onase anomaly jus west of une 
Wasatch Front is most pronounced at periods 30° and 45 
minutes in the substorm maps (figure 4.4) and at 45 
mimubes Por) the storm (figure 4 70)) Au the shorves: 
perlods the phase relavionship would be impossible for 
Pei oem Dysexvernale a 11Lelds., Olle Wher is. cemenpe red 
Chav anomalous internal 4 Vs 100dUCcCed DY VY Une Chase etoma tires 
Arewsecn vo. be COnSisvent with local linear currents 21ong 
tine Wasatch Front. “These linear currents have already been 
proposed to account for the amplitude anomalies in whe zZ 
and Y components across the Wasatch Front. Similar Z 
phase anomalies can be seen near CUB in figures 4.4 and 


4.8 for the substorm and storm fields respectively; this 
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anomaly is located just west of the Southern Rockies. 
These Z phase anomalies can be explained in the manner 

of the Z phase anomalies across the Wasatch Front as 
representing the effect of anomalous internal Z fields 
associated with linear currents flowing along @2 ridge 
under “Che “Southern Rockies, and induced by the normal ty 
frelds Further discussion Of phase taps will be deferred 


to the section covering separated phase maps. 


4.4 Separated fields in the period domain 


As was stated before, two variation events were 
recorded by our 1967 array in the western United States. 
Miewieveids soeuhe subsetvorm lot jsepuecmber Gihaye been 
Separaved sat) LOUr times, Mivuhne ~ame domaine andeay [=> 0 
and 60 minutes in the period domain by Oldenburg (1969), 
andsat T=50 and 69 minutes. by Poratm (Porath, Oldenburg 
and «Gough 97.0). 

The vauvhnor has Separaved ae SOG i soGhs ent, 
September 20-21, 1967 at the three periods 89, 150 and 
256 minutes. As has been explained, the storm was 
poorly recorded in the north-east corner of the array 
ang) ro. Ls in the western Half or the array that. the 
resules of the separation aré most meaningrul. ihe 
eastern stations give adequate control of vhe fringe 


fvelawor che Wasatven Front. “Unfortunavely the Basvelront 
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is close enough to the Wasatch Front that the currents 
under the Southern Rockies contribute considerably to 

Wine Minvernais [held at ine Wasatch Front (Poratner al. 
OT O)MsOstsat ene modelling of the Wasatch Front from 
SvOrm Seperavaon results is necessarily of a qualitative 
Navure. However, Quantitative depth estimates for a2 dine 
Current model of the Wasavech Front apree with those of 
OlGenpure (969) and Porath et al. ©1970). Whe aercement 
Of whe wvresuitS trom Che sivorm wien Those fron the csubs.orm 
means thay tine models pit forward Porm eonductave =< truciemres 
Lieve? Uower mantle sare nov derived Tron a single events, 
ire Workton —Reuiyzel ter al. (197.0), Poracth et a alo 7 Oe 
Oldenburg (1969) and the present author has established 
Uhe wpresence Ob two-dimensivona conductive srructures 
StrLigine @pproximevely north-south under vuhe j0uthern 
Rockies id Wasaven Bron, and Mas placed (MimiGs on tie 


depres, Gimensions and conductivities of these suructures. 


We> wrocbecaved SsLorm tields 
Normal faelds 


lteis now (secer one. Oo) that fields wrth esca le 
lengths greater than the dimensions of the array cannov 
pe separaved. To pemove the inseparable field a “normal” 
hovaZonval tacld ae SaaS cUmed ie sigse AOhevie shigheye Neysy elelS 


vector sum of the normal external field Hee and “2 norma 
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intervie d fel ora due GG Cirrenvs flemane 2 genun an 
vhe highly conducting mantle. The removal of om leaves 
an anomalous internal field, Heo LO Re Seperavea trom 
an inhomogeneous residual external field ee Porath 
et al. (1970) attempted a removal of the normal field 
by fitting ta- polynomial in’x and y by east somares to 


the Horizontal Components X and Y wium the boundary 


Condi tlon 
ax _ oY 
ay aX 


iiiseatcvenply Was UNsuccessnul Since the matrix for calcu— 
ira ol ne. the coefficients became all—condi@t1oned. tor 
Surcieces:; Olwhighner Order than a plane. — Wiem possible 

Ha would be removed Derore  separavion.. However: since won. y 
Splane apocoxmmarron Gan De made to une exvernal tield. 
the f2elds muse De Separaved “first. vO sce 1f tnerp lane 
approximatvaon is approoriave. Une complicated external 
Warus Of Une separaved sborm field for periods 695 150 
and 256 minutes coulda not be represenved by a Constant or 
even by @ plane. Whererore tne mormal Tield could nov we 
removed, nor could the anomalous field be subsequently 
Heolatea and normalized with respect Co the normal Prield 
4m the manner of Porath et al. (1970). It is a4 priori 
eee Go find the mMormal ticld of a svcomm much more 


complicated than that of a substorm. 
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source field 


ine locavion of the current sysvem (sec sece.on 
1.6) which produces a polar magnetic substorm is depen- 
denvwson the Tocatton of the eatroral oval. During amvense 
storms tne auroral oval isi known, to ae equatorward 
(Akasofu 1966). Figure 4.11 represents the relationship 
of the PCuLiWard movemeny Of Tine sauroral oval to whe svee 
for the svorm. The storm or Seprembper 20-21) 1967 has ca 
maximum Dst of 80 gammas which means that during the event 
Cone auroral oval probably moved! scuthward to about 56°N 
lacLvude.) he current, system would them mes lecabedtat 
Soles CLMe aur ine Ene Veil Wal vai 14° of the northernmost 
line of the 1967 array. Figure 4.12 (Walker 1964) shows 
Cnau. wr ne Current 23) Very Wider tiem stawioms, Locared ay 
the latitudes of the 1967 array may be on the southward- 
inereacsine side Of whe fH componenuy curve. 

Wie *k=a1ndiees Of Une Storm tnd care Chav ivomieid 
Was comprised mainly of substorm activity ay Une times jor 
recording. However, wie Source: currents for he eeonag— 
netic storm of September 20-21, 1967 seem to be very 
complex and a simple model for these currents simply does 
not exist as it did for the substorm of September 1, 1967 


(Oldenburg 1969). 
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Figure 4.11 POuUvMerniy Limits iOh Cicer ora lear e 
| in relation to amplitude of D., (H) 
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Ficure 4212 Variepion Of fel onponents Wata 
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ocale length of the external field 


As has been shown in Chapter 1, induction phenomena 
depend critically on the scale length (1 = 27/k). The 
scale length will be estimated from the eperna Is fre ice. 
HO. tie (periods 69; 1507and 256 minutes. 

ue te detined @s the votval horizontva » exvernal 
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Hie parameter ko is Given (Schmucker 1964) py 
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of F,. The values estimated from figures Weis and 4.15 are: 


Period (mins. ) 89 150 256 
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me TOm Comme) flee Cncu 0.65 
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HS 5010 ce eee ©, OOO an 
WwitiMeam average Value Ol 4 01 ~t2,0005 km. 


Thus 1¢ ean be seem that the scale Mengsth or che 


external field is much larger than the dimensions of the 
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Ficure 4a 13 Depara ved. Mielid Seonpeneatis. Or ride 
Svorm aluong Unies Beast -We. Wb Oldlec. 


= 60 minutes. 
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Figure 4.14 peparaved field “compcnents or che 
SUOrm alone three Cacst-wes lm Drones. 
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Pigure 4,15 Separated field components for the 
StOrm alone three Cast—westl prohiles. 
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Figure 4.16 Separated field components for the 
Subs torm elone chres casu—west 
profiles. © = 39 minuves (after 


Porath, Oldenburg and Gough 1970). 
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array. The values are an order of magnitude Lar eer 
than the values of 600 to 800 km found by Caner et al. 


(1967) for his magnetic variations but compare well with 


eine 121m 6s 
5,000 Kin << IF, SOO km 


tound by “Oldenburg (1969) for the external field of the 
SupstoOrm 1 Vseovemnber 1) O67) “Caner etraie (1067) 
Weres Of COUrse, USing Unseparave d fields amd their scale 


lenguns are doubtless relaved to the imternal fields. 


Pxe@ernal svorm tields 


Figures 4.13 through 4.15 Show how the external 
fteld OL The SvOrin Varles On) biree ines acres onune mata y 
aoeperiods 09>. 150 and 256 minutes. © Lt iso, louse unas 
BiewexXLerinal. [helas MUS De ue GO Complex IcUrreiia se vems 
and canhoy be approxamaved by a plane surtace. 

tie exverial Mortiwarcd DOrizenvals tLe lds Ceshows sa 
decrease southward for periods 89 and 256 minutes. But 
Surprisingly the KX Component increases southward for 
Deriod 150 minutes (Tisure 4.17). Inthe paper py ekeiuze! 
eres lo welo 70.) vhe x component, T1501 ot nemo ua. 
field was seen to increase and this was tentatively taken 
to imply a rise in the upper mantle isotherms from north 
40 south. , The expected reversal of the 4 component was 


never seen but this was attributed to lack Of Z Observations 
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Figure 4.18 Separated vertical and horizontal 
eastward svorm 1 reldswau 1 =150 


minutes. 
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in the expected aréa of reversal.” The argument was based 
ons the assumption that the storm external X component 
must decrease southward as did that of the Ssubstorm. “This 
decssnotmholadyfor T=150 mins. feas it does for ’l=39"Mmine. 
andmi=256imins.” It seems possibte thac the euroral soya! 
came very far south and the three-dimensional current 
system expanded so that our array was under the southward 
imereasine i Componeny region of the curvesc! fieure 442. 
The separation was probably not complete since 
Serweulre Can oesecen in Z and as Over the Hast Pron 
anomaly in fissure 4.13. “This may be due to lack -of 
In formavaonvabouu the Tields in this “area. “Ua veeneral 
UUme sex vernal Components ay all tires periods are ssmoo uy 
Veryine oMer the Wasaveh = Mront, “where Une operacion vor vune 
areeayiwas More tt telene. 
The external phases essentially vary smoounly across 
the array, even though at the eastern edge 
(Mere again appear spurious gradients which can be avori— 


Dicee LO extrapolation of Che Maps OulsSide> the “array. 


Internal fields 


Comparison of the internal cosine and Se ire tie Orn 
profiles of the storm and substorm can be made from figures 
Wa? ang §.16. The figures represent Une spatial varciavion 
of the appropriate component from west to east across 


Central lines, in the array. 
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In comparing the two events it must be remembered 
that the signs and magnitudes of the sine and cosine 
Gociiicients depend Upon the starting time of the Time 
series and features of the time series. transform may be 
seeniin either the sine or cosine Berroa Mer Or a th 
Bopha Sinspection of figures Hits, eg. 41> end Yeo 
shows that the sine transform profile of the substorm is 
Very similar. to the cosine transform profile of the storm 
at T=89 minutes, while the cosine profile of the substorm 
Deets IC olewresemDlance To the sine oroi1le of bie. svorn. 
Phe sine profile tor the sverm 1s essentially featureiess 
at T=89 minutes which is most likely due to lack of energy 
iisune sine ctranstorm. 

The storm sine profiles for T=89 min. like the 
SUuCs vor Droid es "show an asyvmmeury over Whe Wace cli mene 
Whiehmcan We sc leary seen aie the z+ and Yy Comocnentsr. 
Poratinen al. ©0700) have shown Gat this asymmetry ais 
fess prominent for Lhe shorver pertods 30 vand) 50 minuces- 
These results suggest that the Wasatch Front anomaly may 
Demduel to 2 ridge Of Conducvings material superimposed mon 
an asymmetric step in the conductive mantle down from the 
Basin and Range province into the Colorado Plateau. This 
asymmetry of 2+ and Yr profiles over the Wasavca Fron. 
persists at periods of 150) and 256 Minutes, sSucegesling 
that the woheaval and step are deep structures and may be 


continuous with the highly-conductive mantle below. 
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Figure 4.15 shows profiles at T=256 minutes. 
Inspection of the sine profiles and the profiles “at 
T=89 minutes (figure 4.16) shows that the Xo, Yo, Zp 


COMmpCNeNnUs are very similar. Both the Xe and Ye Components 
are very smooth and can be approximated by a curl-free 
plane. Both 2. components show little east-west gradient 
which suggests that ie may be Set equal to ertnen (Poratn 

ev al. 19/0). The procedure to obtain be and Ys, was 
then carried Ouv ion the sine transform components av 

T=256 minutes in the hope of amplifying and focusing the 
Wasetvch Pron anomaly, ‘Inhesresults Of Unis analysis scan 
Deseeen im the maps of figure 4.70. “The maps ot Zo shows 
tie Wasavcherrony learly bul ay Seems TO We woroauened 

ONE COsUune Ppreaver eCxventl Of tthe anducedscurrentssunder 
ener Colorado Plateau av. Giis lone period.) No tstave iis 
another way, the Wasatch Front anomaly at 256 minutes 
period would De more intluenced by the Base Brony anomaly 
tien savas UMM EIUU es. es Lerma pos Lia shows large 
Pradvenus roughly outlining vie Colorado Plaveau. |Far— 
Cilcular Unverest resiges in the larce gradients an Zea ala 
the western and southern areas of the array.” These 
Bradientvs may reflect the Doundary Or vnhe Colorado Plateau 


and the Basin and Range province to the south and west of 


the Plateau. 
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Me ceneral patvern, of the Ls PiLeld. ie pie 
mur cyois ther array imereasing Go. the, south, supporis 
the idea of a north-south isotherm under the Colorado 
PisGeeuerasinge Up inte, the Basin and hanee. FO we le. Soul. 
iiese results, seem Fo support the aie Chie lar sine ace 
eviain, (1970) mentioned earlier, 

The Y; sine maps of figure 4.19 show the Wasatch 
Front anomaly very clearly and the trend of the gradients 
MS hoceaved west Of the trend of the z+ CeadLen sme wale 


Seely aoreeme nin Wl oly walle resultseOn Pore um em, aise Olio). 


Comparison of the ce and: VY ~sine mans at 256 manutes 


Ai 
Giteure 4519) shows that the Wasatch Bront. anomaly 1s 
Peeankyecemhianced: Diy jaie remova h.OL che normal  aanvernat 
feds, 

ihe spectrum of the.Z component «of the ssvorm 
Gircure: 3.5 ) reveals play most Ol Lhe power minetne ss corm 
is Concentrated near T=150 minutes and 110 minutes - 
Kigure 325 salso shows vne relative -amplivudes ol the 
three eomponents H,D,Z4 at various. periods. Tt seems that 
Ghe ratio of the amplitude 27D and 2/7 is greatest aise 
Se 150) manuves, Therefore 150 minutes. is anjideal period 
at which to separate the storm because: 
(a) there is sufficient energy available in the Z component; 
(blotGhe ratios Z/b and 2/7 are favourable- 


The profiles of the cosine and sine Z+ and YS com- 


ponents at 150 minutes (figure HA) again show the 
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asymmetric fields over the Wasatch Front anomaly, but 
much more clearly at this more favourable period then 
at 89 or 256 minutes. Maps of the Z> and Y, amplitudes 
Giese 400) eive a ecod representation of the Wasatch 
Front anomaly, while little appears in the sparsely 
Observed= caster lial (of the array. Tne J0call on ot, 
the maximum of Yr iS again Seen to be wes Ore tie Ze 
Maximum over the Wasateh Nront; this supports the location 
of an anomalous Line Curren minder yo Neu Wasa oem i ieOmin mm le 
Xp, Ze Piase Maps seem to be yery smooun and) Uneven tiuk. 
while Xv> 2+ phases show definite abrupt 
Cranvecs iil sven. IMS peer en (Ol auine Xs Ze phase maps shows 
WieiwewMese solace ce langes rong ouG line Theewes ver and 
southern boundaries of the Basin and Range wand the Colcrade 
Piscedu.. This Outline CcoOlncldes Witn tie soul tier ooserved 
eater Ors Lie Zs sine tramsfLorm maps, eat T=256 minutes. 
In summary, the Wasatch Front anomaly shows in maps 
Of the internal parts of Une svormorteld components. spain 
ticularly at period 150 minutes at which the storm spectrum 
has vebunadant energy. “Whe resulvs coniirm whe presence or 


an internal conductive structure in the position and azimuth 


imdicated by analysis of the Separated subsvorm gields. 
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4.6 Conductive structures 


Two types of conductive structures can account 
for most local anomalies observed in Geomagnetic Depin 
pounding investigation (Sectiom 1.5). 

(1) Surface anomalies caused by high conductivities 

iM erus val rocks lose) fo the surtace, Im many cases 
wmelated Vo -decp Dasdns of morcus Sediments. The anomalous 
internal field should be out of phase by 60° to 90° with 
the inducing field since Ohmve resistance: is more important 
Pia aucudnre wire acvaricer, 

CO; Undulations of “tseocons"” “surtaces Of equal, conauc— 
Pivity) 0 a Conducting Medium oF darce exten.) lhis 

Pye OmnoOe LV Lee represeiicauinvic mol ce DmCONGUCiI CS 
Sucvletuces ia the Harv. Theses undue apions dave  conduc— 
Dive Ss NG sOtmenstOns,.Characver sual Oh wUppei amaitte 
SULUCGTUCe, Such thav (IS6li—inductance conurols ule anouced 
Currents rather than resastance,  Areuments presenvedein 
Chapter 1 State thau the phase ditterences berween mormal 
and ee wou fields associaved) With such Structures are 
small (less than 40° typically). 

Porath et al. (1970) estimated ihe rea leconducuiyiw, 
for che sem -circuUlar Upneavel Under Vie sOouvnern hockies 
from phase differences between the normal and anomalous 
fields. This upheaval has the same response as an isolaved 


cylinder (Kertz 1960, Rikitake and Whitham TOGO) SW i ath 
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associated induction parameter c = r/4nwo which can be 
estimated from the relative magnitudes of the in-phase 
and out-of-phase fields. The radius r of the ey Lander 
can be obtained from the maximum amplitude of the 
anomalous horizontal field W = Re (Sechmucker 1959) 


irom the expression 
P= 2yv W/E 


with ft tne fraction of the dipole moment of a perfectly 
conducting cylinder, which can be estimated from the phase 
differences. Thus for the Southern Rockies z = 360 km, 

W = 0.175 and f = 0.6 yielding r = 200 km. For the phase 
angle BOs Ce eo Cr Uy Wil tel C eves eee OTC eva tmGeemcas 
OS ide CaM. ©Or tne Trdee Under= tile DOUucnermenockiles 
(Peorath ¢«¢ aie 197.0).. 

(here Cxisy Several reasons Why Une vwoO Morin — 
SOuUchmanomalies. Can be imverpreted sas Deine ass0ctaved 
wre laree, deep conductors: 

(1) Whe general correspondence between the anomalies and 
heau flow resulcs (Retezel et al.) 1070). 

@eihe tneability of sucrace anomalies. TOmaccoum Or 
Diace Cirterences, Of abou 0" between the normal and 
anomalous fields (Porath et al. 1970). 

(2) The persistence of the anomaly ab very long periods; 
this work showed Dereisvence Up vo 43, hours, Ba Reitzel 


etal, (1970) showed persistence of tne anomalies in the 
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Figure 4.20 
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eq daily variations, which have Negivzi bie energy an 
Harmonies Hieher than the fourtm (rT = 6 hours. 

The anomalous fields produced by lateral inhomo- 
geneitvies in the conducting upper mantle are a second 
order effect superimposed upon a much larger normal field. 
Me laveral conductivity inhomogeneities are often inter— 
preved as Upheavals in the mantle. Li these upheavals 
are modelled as an isolated cylinder then neither the 
line current nor a two-dimensional dipole resulting from 
TAS I NdUCT TOM are physvcadiy wealisipae Tor two Mmeasone. 
(1) The line currents allow for no return current. 

(2) The dipole postulates the existence of two anti- 


Davee l CUrrenys In vie Woheava ly. 


Rie eure remy, COneen urs biol One Gileme We VeaiGe Gao mimi 
Of wyene upheaval can be approxima ved by a lane icurrenn, 
Wolo view GitutwoOn Or Che eGurerenD CensiuyOllee voles tae 
Ofeuhe Upheaval Ts kquivalenuy to a return current under i.he 
tine current. Thus 2 semi-cylindrical upheaval or aniinice 
Conductivivy effectively has an’ image of the forward eurreny 
under the upheaval (Rikitake and Whitham 1964), so that 
the half-cylinder has the same induction response as 4 
conducting cylinder embedded in an insulator and produces 
PreEpield of a tworcimensionad cdipose (fipuren4 20)” 

The field of a semi-cylinder on the surface of a 


perfect conductor can therefore be approximated by a tine 
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of dipoles at the depth of the undisturbed surface of 
the conducting half-space or by a line current near 

the top of the upheaval. The line current then cives 
a Maximum depth to the top of the wpheaval. The two- 
Gimensional dipole model is used to estimate CeO cine 


the undisturbed conductive surface near the ridge. 
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where a is the dipole moment of the two-dimensional 
Gipole. MM and N are the in-phase and out-of-phase 
components of the anomalous field respectively. 
Equations 4.3 give the depth relations of the hori- 


zontal and vertical fields Cespeccoryely: 
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BV 52) Ours 
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Depth estimates were obtained from the spatial 
half-widths of maxima of the separated magnetic fields 
alone three centrally located east-west profiles, Figure 
Wea egives vie amplivude protiles: or Yr and z+ ear. Se 
menuces. — Une depvuhn esvimnaces were made only tor wcne 
Wesel Prone Ns GSivaAved Drevlols Ly Folly une ne came 
transtormeaay | = 250 minutes Was wt appropriave to, teolabe 
the anomalous field in the manner of Oldenburg (1969). 

Al sothner depth estimeacves were made from Uniormala zed 
Separaced magnetic variations, Comparison of the asolaved 
Go the non-isolated fields Shows that the values are wirhin 
thee 2O ter rOC® acSOClaLed Wetm Cac ectamalens solnce. ule 
depth estimates were made from unnormalized fields a base- 
line representing some portion of the normal field was 


G@tawn. Depths to a line current and, intinive cylinder 


are given in table 4.1; all estimates are taken from 
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magnetic anomaly profiles like tigure “22 over stne 
Wasatch Front. Comparison of the results above with 
those Of Porath et al. (1970): can be made by inspection 
of the depth estimates they computed for the substorm, 
quoteds in. table 4.2. 

Timevdepun estimates ©f the postulated scurrents 


are consistent for all profiles and interestingly enough 


fom the very Wong period also, elving the vaverave maximum 
depth to the top of the conductive ridge under the 


Wasatch Front as 
LV km + 207 


Ol@enbure (1909): computed @ Similar table to 441. for line 
Current. Ceprn estimates Por tne SUSTOImM, ant ound or 
average maximum depths to Wasatch Front 110 km +t 20% 
Resy Prony POUR 20y, 

He also reported that depun esvimates siren similar 
prefiles on maps of the sine and cosine transtorms som 
Zen and Nes and om maos goo Bee) lew and Re la 
also @ave maximum depihs of the lime current waich) were 
Within the Limits of the unmormalized depth estimaves. 

The results of Oldenburg (1969) and of the author 
show that good estimates to the top of an Upheaval modelled 


as a line current can be made from the unnormalized separa-— 


ted eleiaey mhe consistency of the depth estimates at 
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Table 4.1 


Depth Estimates of Wasatch Front 


Period (mins. ) Depunt( from half-width )km 
Line Current Cy lander 
89 110 Oy 
150 ah alay/ 250) 
256 V4 2s 
alee a2 


DEOU Rn esuamales waEOm sUbSUCrM er Velde 


Depth(from half-width) km 


a eae Dane Curren’ Cy tanger 
Southern o2e5. OSE 180 370 
Rockies BOs 0.9 190 4OO 
| 89.0 Omg 200 410 
Wasatch SOS Lgl) ALEO 250 
front Beh J0) 2 5 260 
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long periods corroborates the COnciIusien of —Reitgzel 
et al. (1970) that the anomalous conductive structures are 


deep and may be regarded as undulations in the conduc- 


ting mantle. 


4.7 Models of conductive structures 


Porath, Oldenburg and Gough (1970) have developed 
a fave—svep method vo obtain models of conductive 
structures that cause two-dimensional magnetic anomalies 
observed across an array. Whe five steps are: 
(i) A line current model is used to estimate maximum 
depchis to. ume TOps of conducting Upheavads. 
(ii) A two-dimensional dipole model is used to estimate 
taekcepuheve whe Undrstirbed Conductave sSuriace, Mear Time 
Pra ece. 
Git) A semi-cylindrical upheaval on the surtace jor a 
Derrectay conducting half-space Vs Uriled, starving warn 
depen and radius as eavimaced from (i) and (1) ey wick 
mecessary adi ustmenvs to improve whe Tic.  Wne frelas 
aeseociaved With such @ model can be computed by the metaod 
of conformal mapping (Schmucker 1964), 
(iv) Asymmetrical anomalies can be fitted by a step in 
the surface of the half-space, and combined in a fivst 


approximation with the halfp-cylinder model 21 mecded, 
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(v) Once the dimensions of the required conducting ridge 
are eStablished phase differences between the normal and 
anomalous fields can be used to estimate the true, finite 
COMCUCTIAV It ys 

Mis vhesis has compleved sueps (1 )) and (406 Whe 
Weiter dad NOL preceed tO Steps (i 11)), Cay) andy because 
the Vesulvs Of tie storm separation CO nou allow cuany.— 
tative modelling. The separated fields over the Wasatch 
Front have been compared with the models developed by 
Porath eval. (1970). As was pointed out preyucusly=the 
anomalous fields muse be mormalized to obtain 200d estimares 
of Une CONndUucCUIVILCTY Using the piase ditterences Werween che 
MOLMawh wand Aton lous sme lOSse. | Ine mess son ne sLOrm 
Sseparaclon Mave shown the Pvelds vo Be Coo Complex vor be 
normalized joy present methods, so that no quantivavive 
egtimave Of the COnductivaty was Obtained from tie Svorm 
fovedasr. 

‘Porath eb al. (1970) model the two anomalies by 
fitting the observed anomalous fields to the response of 
semi-cylindrical upheavals and steps in the Surrace ond 
perfectly conducting half-space. The responses were 
computed by conformal mapping techniques (Schmucker 1964). 
Figure 4.23 shows the steps involved in approximating the 


substorm vertical field anomaly by a conductivity model 
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Fagure 4.23 Model aT Veted Go norma lavzed snome lols 
Verbtucal Tileld vor supe vomme: 1 =e 9 
minutes lafter Porath, Oldenburg and 


Gougw 1970). 
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Figure 4.24 Model fitted to normalized anomalous 
eastward MNorizontal facld vf subscorm 
ac = 00 mniges (etecr koOraun, 


Oldenburg and Gough 1970). 
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WarGer Perauh eb al. 1070), Figure 4.234 16 4 Semi 
cylinder on an infinite half-space whose dimensions were 
derived from results of the line current model and the 
Beene model. Figure 4.236 shows the response TO 7a 

step model compared to the observed anomaly minus the 
field due to the model in figure 4.23a observed over the 
Wasatch Front. Figure 4.23c shows the result of super- 
position of the step response to the cylinder response 
model, which then gives a good approximate conductivity 
model to the observed Wasatch Front anomaly and the East 
Front anomaly. This procedure depends heavily on getting 
a 2000 estiamave of the Haste Wronusconducvivi py ss VeuecuKe, 
Because Such an €Stimate was mov obvained by the writer 
hor une anomalous Storm faiclds this medetiling procedure 
COuULaGsNoOc bP -compleved. The work On une sseparaced: svorm 
fields’ terminates with 42 qualitavive fic of Che unnormai zed 


Ye onde Zo anpli bude COmponenvs Or Che storm au tT > OU inutes 


i i 
as profiles over the Wasaveh Prony. Figure 4.25 shows 
Evo ce wo anwne  limins On Une WepUm esuimaves Yr and Z+ 
Pic une theorevical responses, of ule YT and zy components 
over the Wasatch Front as seen in figures 4.23 and 4.24. 

Aithough the lack of information from the eastern 
half. of the array and complexity of the external field of 


the storm prevent completion of a quantitative interpre- 


tation of structure from this variation event, the 
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qualitative and semi-quantitative results secured for 
ohne Wasatch Wront Structure are important’ in giving 
independent support, [rom 4 second inpus field, ua vine 


conductive structure deduced from suudy of ene ssupsvorm 
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CHAPTER 5 


GEOPHYSICS AND UPPER MANTLE STRUCTURE IN THE 


WESTERN UNITED STATES 


Geomagnetic deep sounding results from the 1967 
array in the western United States (Reitzel, Gough, 
Porath and Anderson 1970, Oldenburg 1969, Porath, 
Oldenburg and Gough 1970, and the present thesis) 
provide tle opportunity for correlation wich obner 
Seophysical results. Such results anclude refraction 
data in seismic P velocity (Herrin wnd Teeearc  oG2. 
Archambeau, Flinn and Lambert 1969), seismic time-term 
studies, (Cleary and Hales 1966) extensive heat flow 
measurements (Roy et al. 1968a, Roy, Blackwell and Birch 
1968b) and conductivity-temperature relationships (Tozer 
1959. Hughes 1955, Switt 1967, Hyndman and Hyndman 196307, 
Warren eo pal. 10695 Goch and Poraga lo 10). witli tae 
hope ©f e€svablishing 42 geophysically consistent areal 
pavternm Of Upper mantle lareral’ anhomogeneities: in) the 


wesvern United Svates. 


5.1 Seismology 


Herrin and Taggart (1962) applied a computational 
procedure to seismic data in the P, range from 89 
stations to determine the epicenters of the Gnome explo- 
sion and the Lake Hebgen earthquake in the United 


deates, The investigavion resulved an une finest 
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clear evidence of significant regional differences in 

i Velocity. in the detailed ahalysis: Of rercvonad 

Var lapions in a Velocity. tie iam phase is regarded as 

a head wave travelling in the upper mantle just below 

Tne Monerovictie CisGoubanul tye G distances less than 

2000 kms implying that whe i VelLOCiIuy= COes Om. vary. 
Systematically with distance from the source. Interval 
velocivtes Were €Ompuved betweem Stations alone radial 
fines from the source by dividine the Gistance between 
stations along the radius by the difference in errival 
times at “the two stations. Whe validity of Cis method 

iS Ose O Onl UINesCONUI VANOn =U ait ue aCieUstronleriin aciciaie s cumin 
une same vunder bDoUn svalions.,. This teonditiom would thus 
oreveny the devection of the Wasatch Prone where Gis 
Pesulus SUSzeESst tna phe crustal wnickhness schenges DPeuweoen 
the Basin and Range Province and the Colorado Plaveau. 

rhe assumption uhat. invervyal yelocilmcs pprosimeve sune 

oe Velocity ay the center of the spread allows the results 
to be contoured. Ricure Sal TiIusvraves wie Compr es carmOia! 
velociuy at the Cop Of tie Mant Me according vos the 
assumptions given above. The attenuation of amplitude 

im the f£iret arrival from whe Gnome explosion Was seen 

to correlate With regions of low a VeUOCisG Vw 1 Vieesa barons 
faethe dept of a low velocity layer in the upper mane 


then may well account for the complex ae VELOCI Dy Pavvern 
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APPARENT Pp, VELOCITIES IN 
WESTERN N AMERICA 
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(Herrin ond Taggert 1962, 
Herrin 1966) 
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in the western United States. Herrin and Taggart (1962) 
maintain that these variations in depth and thickness 

On Ehei velocity nimimum are probabdy: Pelated Go 
Ginterencessin terrestial” heaturlow and thus to pe 
thermal, pradrencs, LWousiie ope present an Ehe voper 
Mantite. The averaging of ae Veloc up ies alone 2 prot ile 
definitely limits resolution, preventing the recognition 
Ore cuchi re lava yely narrow features as the Wasatch Front 
and (southern Rockies Upper mantle ridges. 

Cleary and Hales (1966) aneadiyZed= the travel Times 
ofvPewaves from 25 earthquakes! to ‘stations an Norch 
America, From sources at angular distances between a 
and TOOL 5 by a@ileast Ssqmares technique similar vo ‘che 
"tame term’ method of refracviion seismology. Deviations 
from the Jerireys—Bulten tables were separaved irom the 
errectseor Staulon sndesource residuals.) ihe dirterence 
from the J-B tables were @alculated at 20 antvervyals and 
Cictced tora smooth Curve.» Ao regional rend was snow 
ial the Stavuuon arrive tames Wilh errivals an basaneand 
Range Provinee being Up to a second lave and vhose= in 
fhe central’ pare of the Unived Svates Up vO a second 
early. 

Cleary and Hales (1966) decided that data should 
be selected initially from stations on the North American 
continent to Keep tne eEfrectus Of Mocatlon errors vo a 


reasonable minimum. Their reasoning follows: two 
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stations at opposite ends of the continent at a distance 
06 from the souree, with a mislocatton of 0.1 produce 
a difference error of 0.7 seconds and a distance of 60° 
with the same conditions producesa difference error of 
about 0.5 seconds. Earthquakes are chosen that 1) were 
shatlow-locus, 2) had 4 first arrival that was sharp 
enough to be recognized and measured precisely at vir- 
tually all stations, and 3) were distributed in such a 
Mander that the data from most statwons would: convain 
information from an adequate sample of distances and 
azimuths. Ihe Stations that provided dava Tor the 
ana ly sas were: 

lt. LRSM Svaviron and Vela-Uniicrm Observatories 

Zw MWsoS oval Lons ny ties Uso. 

ie SU cuOuLSs NetWork Stacilon 

4, Caltech Network Station 

bn Canedtan soualvions 
AMore cis tom Ot 01 Sec. was attained for le Wire corda, 
A version of Bolt's ID1 computver program was used to 
@aleulate the epicenter, or give time and focal depri 
for each event alone with the residuals from tne iB 
times. The Nortm American station residuals were 
extracted for analysis from the above residuals. 

Me method of analysis as based on evaluarving 
Give equation 
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at each station r and earthquake s, where oe Ls oie 


observed residual from the J-B tables: Ang? assumed to 
Dea Tunctlon of distance only, ws the everaze “dit fe— 
rence from the tables at the distance of Station © from 
earvaqueke is; bo, 1s) Ue s bac Ofueces dd alunite ase. 
perturbation in the travel time induced by conditions 
iene Mer siOrncod oF station sms sand qd. 1 a2 perturba— 
Peony Caused by Condivions peciiiiver tO eartaquake Gieaas 
IS! CWC 25S. OOMS Gens Gl elena ie IEMA NS NOONRete ere! leis telete 
data, As Menvioned above the analysis 4s samiler co 'eime 
term" with the difference residing in the parameter d.. 
Pitceaks NOv really a cOurce  Geicn OU ele OU lGmine hes pieOmeily 
DewCeal leo an varoitGraryebese Line COmrectmonemade seca sooty 
Decause Only a) SseclecCled Tou On Suellona, Were wits ec 
thesmean sol the residuals of thas 2roup May crit er irom 
Zew er 

Tre results Meom che: analysis of “sna low-nocus 
earthquake data are shown as station P wave residuals 
once Hales 1966) over the Unived States in figure 
Bee, There is indeed a marked regional urend wih 
large positive residuals in the Basin- and Range Province 
and large negative residuals ian the central Unived staves 
with apparent extensions of this negative residual trend 
north into the Canadian Shield. Figure 5.3 was vaken 


froma Avaeere by Doyle and Hales (1967) where a similar 
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TPIS OOS Sig 2 SUCH UO Al igi eiOmUlewibiske jail: ie ehiahek yey ll mallets 
over the United States (Cleary and 


Halles 1966) 


“1.9 


LRSM stotions 


WWSS stations 


Coltech stotions 
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Li ge 


Ia benbhigkes Baye.) SbavLOn resmdUalsS an omer riy a ames 
over the United States (Doyle and 


Hales 1967). 


APPARENT Pp, VELOCITIES IN 
WESTERN N AMERICA 

Adopted from 
(Herrin ond Toggert 1962, 
Herrin 1966) 


def A. 


analysis was applied to S waves as was done by Cleary 
and Hales (1966) for P waves. 

Cleary and Hales (1966) based their analysis on 
an earth model that had a laterally homogeneous lower 
mantle overlain by a laterally inhomogeneous lipper 
mantle and crust. They have also assumed that: 

1) the average travel time of a P wave is con- 
trolled by the homogeneous material through 
which the wave travels for the ereater part 
Oty Pibiotce weyehngels 

2) only inhomogeneities in the upper layer 
produce variations Prom This averece tame 
end are {constant Wlthian Temp rescri ted 
GistancevwRawse hor a spect lier sire vlonm ra AG 
seurce. Cleary and Hales (1966) state that 
CNe COnsisSvency Of Vier eesti Ts suppers 
Giese assump eLVOnsS Or ene Tres morgen, 

Thesuse Of en errey Of Seusmoneters On a wesc tonal warier 
than a continental basis in a Cleary and Hales "time 
term” method of earthquake analysis Offers promise of 
providing a seismic technique capable of comparable 
Wateral pesolution of features to Ge array studies. 

Archambeau, Flinn and Lambert (1969) attempted 
to determine the upper mantle structures that underlie 
particular Peneaneat al orovinces through the use of 


modern computer-based techniques of processing and 
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Figure 5.4 Observational profiles used by 
Archambeau, Milann and Lambere 
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Make nee Ss eissniice ace co Syetvematically isolate and 
identify mantle body-wave phases. they Sbave nap 
they have been able to obtain 2 large amount of relia— 
bie iantormapvon from a relatively small number of 
observations along a single profile or linear acreaye, 
becawse their Techniques allow data to be obtained 
fron Maver arrivals. uhus Giving complleve travel time 
andweamplivude Gusvance Curves, Pariwer than jus firer 
arrival branches. These authors obtained compressional 
velocuty models for the North American Comtinent watch 
distinguish Gdistinet mantle structures for two tectonic 
regions: (i) a basin and range or rift zoné mantle 
Seructure; (11) 42 platean Mountain structure whicm may 
SOW Urner subdivision as more idave ec onmerayas tobe. 
Archambeau et al. (1969) used the compressional 
pases mrom une Shoal and Si loy nuclear explostous a1 
Nevada and from the Fallon earthquake (which occurred 
at the Shoal site) observed along HOUT Oe suCa Lys iae ad 
profiles across thie DAyoLesrapine Drovences, Of erige 
Conminental Unuted Staues sll eure oni) elon cs pitot eles 
were expected by Archambeau et al. to provide a good 
sampling of the Basin and Range upper mantle, the 
Colorado Plateau and Southern Rocky Mountains, the 
Snake River plains, the Norvtmern Rocky Mountains, and 
the Sagem Basin and Range provinces. Because of the 


limited range of these stations and the probable tran- 
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Pare uieS 5. Velocivy {pret tics ound Oya renatbeaun 


Plinn and amber: (1969) in four regions 


of the United states. 


Compressionol Velocity, km/sec 
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sitions in mantle structure alone these profiles 
Archambeau et al. (1969) only hoped to observe lateral 
variations in mantle structure above depths of 300-400 
km. 

The salient results of this etudy ean be seen 
insti gure So, pervinent To this: d1ceuce on an the o 
Zones which Jie in the depth range in which lateral 
variations in the upper mantle properties occur. Major 
changes in velocity can be seen at 150 km, 400 km and 
650 km. Archambeau et al. INPerprevea. Wopem manu le 
SUrUCTURe appropriave GO Une Souvhern Rocky Mountains 
amd Colorado, Plateau provinces, as cCOnMsieting on 2) Low 
VeElOCIUYs ZONe iCa0ped pyaa nveh VeOlociiny Mad Zener = Wi on 
Varia cons in oun the Lid Zone and owe velocity ezones 
from province VO_ DEOV INCE ine asin and: hanwem tl cmsa 1a 
VOlnaVve sagvery Una dvd sor none at. ani wall aa binorma lly 
Vow Velocities exvendinge from or near tie pase sol the 
Creuse to 150 km. A pronounced low-@) Zone 15 Said uo 
correspond to the average low velocity zone depth range 
for the velocity models. A major resulv of Archambeau 
et al. (1969) is the existence of laterally variable 
Pome avor the — velocity i tle Upper Mantle, wath rapad 
changes at both the top and bottom of the zone, as well 
per tne decrease of low velocity values im between. ini 
result suggests that partial melting of Mantle maveriak 


is occuring in this’ depun range. 
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Figure 5.6 


NAG) 


EVMOMUOnWOR GoReniie sa anGs at nem feleci ries i 
COnCUCUAyivies, aiver Hyndman and i yiicmian 


(1968). 
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This method of seismic data analysis along 4 
profile still dees not have the capability of resolu— 
tion necessary to accurately map the boundaries of 
the upper mantle provinces, because resolution is 
limited by 1) the averaging process of the Fe vyelocacy 
resulvs and 2) the spacing of ‘the stations on the 


prota les 


Deve Conductivity and, Thermal suruccure aso, tunction 


of Time 


Hyndman and Hyndman (1968) suggest that most of 
thewelectracal sconduclivity anomalees Opservedsam tic 
North American Cordit leray and partaculerlyethat sos 
Vaners eb al. (190/), whieh indicates chat tien conduc— 
DiVveney le wlareely sre stir CLed) GO a. col kmevbave rilavaras ce pul 
of 20 to #0 km, are related to highly conducting material 
associaved With the existence Of waver—satvuravederoeks 
im che Lower crust. Hyadman'’s: arcumenv is) based orm che 
development of a geosyclinal crust. This development 
Consists of three stages representing ditfereny 2nstants 
of geological time. The results of this development 
with respect. oO conductivities. and. temperatures <at 
depth are shown in figure bwow, = lhevalternavive, posei— 
bilities of electrical conductivity anomalies being due 
to lateral variations of conductivity in the upper 


mantle or being due to water-saturated rocks are 


ee a 
LP lp igamre: sfegiane web: 
sade to Gottedease! ent, a 
Tey) wel antes etd qat €F 

af notzwlet> , -epansce) 42 
efteule? % wie’), saeeray 


eis of waehdae 6h td @rieagneaial | 


. jart@ves (SOO) Wei? ey Pact 
ee tan tutgete odd | 4 
Weette  e Leo? nasties Ga “ 
bier: roof) atdngata ete oeyerey Vaeogs none _ 
degeee 4 wok Ci 6 (pith qtayred wl valved 
be | 3 poem 4 basatiws, agg yom, OF 02 0S : 
eae" hottie todew. 4g (6 aie eeeae one shad tiny hee 
wh. ne 1S § five gog «Peeing tame “aes ati 
(itinutayas <tet jan inalP ayaa Ms aye qazs 
ears bert: rues @riiirger * ORSase yeasts a 
gil OWOd 02 ete Salads sf" ..8 pe 
Je meni (ee ant ee 


wfatinwg 4) tacriattz aot re sth ma 
eet guied -w4i Lenoge, ye i 
7 eae Tilinky ne dot 


: 


nog 4 odd. Ae OVE 

< 

,* ul a, 
| a eae ae De 


176 


dependent on the age of the area Deing Btudied, Thus 
the interpretation would depend on the region being 

in the "Young Geosyneline", "Intermediate" or "Stable 
paield” periods. The importance of the age of tae 
rocks is caused by the decrease in time of the effects 
of water saturation on electrical conductivity. 

Gough and Porath (1970) suggest a common origin 
to une radge producing a heat flow anomaly and to the 
BIGSe Producing a magnetic anduction anomaly. across 
une SOutnermn Nocktves: (ReitZeleer sal, (070m -oracinc mal’ 
LO7 0). ii “Une Origins are Commom “them tThewridces, muse 
bewom COneidéerabile, aze. ae.) How tiermal sconduCcLuvaurcs 
LiauMe Crus. and! Mantle Pilye wa Time Consvont vom melemet 
Hi ormey LOW Ort cinapine eal) depiul a wine “Order sot 10° 
years. Gougi-and Porath (1970) comeluded then that if 
Tie adOve ASSsUuUMpULONS Were true Tien Tne surucuures 
could not be less than 2 x 10° years Old “Uhas = conciu— 
sion rests on the assumption that the induction anomaly 
is caused Dy suructUre, probably thermal; “ane vbe upper 
Meine Le-. 

The writer suggests that the two main north-south 
anomalies Wasatch Front and East Front cannot be inter- 


preted in terms of water-saturated rocks. 
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Deo) Lee Pica lL Conductivity, Temperature and Heat Flow 


The “electrical conduction im the neaer-surrace 
Fock formation is mainly SlocLeOlyiie 25. a) ecm vot 


saline solutions filling pores and cracks. The uncon- 


sOQlidated sediments have oooh siepen e between 0.42 and 
1 (ohm ay in contrast to dense and therefore poorly 
Coudvectin= Aeneous Tock (0. 0l=0.00N comm a). The 
Grust #3 essenvially 42 poor sconductor Since rocks beeome 
More wan more insulating Under ocescire.) a lew newconcuc— 
tive Dy (OCs so1Se (Stein ar ereaver depths anc enlist ise 
ie argely avtribuvable to Amerease agnstemperacure as 
depvuhe 

Ringwood (1966) produced a "pyrolite" petrolo- 
fuca model as al A0prOkiMay lon WO. vie eConle Osa tr VOLmone 
GS WeOeIe iileyenedke > “Inne: siMocke iy Sehy assis. exea snare. 
Gewteliy aaa chemical condiyions. yield Ve slass a 
chemical composition 75% perodotite, 80% olivine, 
(Mg ,Fe),Si 04) 207 enstatipe Mewon 02) SnOm25 (ebase ioe 

“Ihe main constituent of Vpyrolive! olivine was 
impestagacea by Hughes (1955) With respece co als  con— 
ductivity a6 2 fUNCT LON Or ~vemperauine and spueseurec 
Meee les3 kk and Pols % sndpetweens0mand sot ecbers 
imewampient pressure at 30 skm depen. 

Hughes' results indicate an exponential law wake ve 
the conductivity of olivine; such a relation is typical 
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where A. is the activation energy for one of the 
Ssemi-conduction processes, IMPUra ty Ie rinas 1c .0r 
jonic. Tozer (1959) has shown that the predominant 
semi-conduction processes in the upper mantle are 
intrinsic and onic. He also states that the A,/kT 
exponent is large and the dependence on temperature is 
Very svrone, While the error™introduced: into the ten— 
perature calculations by thespresence of tovher phases 
amounts to only a few tens of degrees for quite large 
"impurity" contents. Therefore the precision of tem- 
peravure deverminay.on from electrical conduc iv iiye4s 
highest in the upper mantle. Tozer (1959) has derived 

temperabuee ay depth trom the conductivity models or 
Eahirie and Price® (1939)* and "MeDonald (l95/)] “The 
eonductivity-vtemperature relationships are svill net 
poo well defined in the upper mantle, mainly pecaucs 
of sneeroainey With recard: (Osvne Lavyawice™ CONTENU ROLe ue 
Ou Valine . 

Tt the surface of the earth and that of the 

the highly-conducting upper mantle are assumed to be 
isothermal surfaces, then the coincidence of heat flow 
anomalies withthose of electromagnetic induction should 


ZL Ve information concerning the conductivity-temperature 
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relationship in the upper mantle. This was the approach 
Vakem by Warren 6b al. (1969) to Conpaere tieir Meat Flow 
measurements along a profile coincident with Schmucker'"s 
(1964, 1969) profile of magnetic observations in Southern 
New Mexico, Arizona and Texas. For the ecineiienee or 
the tnermal and electromagnetic ere tO Nave meaning 
in regard to temperature, each must be shown Leaeoendenuly 
to arise from the lateral inhomogeneities in the upper 
Mange.  sehmucker (1904) has @iven @ reasonable inter— 
pretation of the Rio Grande anomaly in terms of lateral 
Conductivity inhomogeneities of the upper mantle, and this 
nas becn Ciscussed an Chapver 3. 

Roy et al. (1968b) haveshown that major changes 
in heat flow taken from 38 bore-holes in plutonic rocks 
aeross Tne seonulmenval United Staves cen woe va 5uripuved 
GO) ehanves ine tne source sot neal av depths wyaciim 
the upper Mantle. Whee analysis vakes sccount jot G16 
Dea peprouluced Dy eradtoacti ve decay in the plutons by 
Us ine the linear heat flow relationship Q = a + bA, 
WieremG de. une meau 1 lx sab wUNe Uciace je is ema lo 
SGU, meabh sprOducuLOnN On es Uirkace plutons, Dune Certs 
to which the radioactivity measured at the surface 1s 
constant and a the heat flow due to sources in the lower 
crust or upper mantle. With information about radioactive 
element concentration at each site of heat flow observa- 
aon ae lloras We for bh can be estimated. | Boy ound what pb 


varied from 7 toll km across the United States thereby 
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allowing a standard plot of @ = a + bA with coordinates 
Q@ and A to be drawn. - The Invercepi Of O62, would ae— 
present the contribution to Q@ from the lower crust or 


upper mantle. Thais shows that heat flow variations 


eoscel. Ona Marge scale ere due mostly tO Variation in 
pie Lower “rusts and Upper mantle, 

TRencombined Tadiweacwuy ity sand, Mea 5) tow mea cure — 
ments at the 38 localities show that the continental 
UniGed  Staress has ieee definive heat Clow provinces: 

1) eastern United States where b=7.5 km and 4a=0.79 

u cal/em*sec CHlPoU2)G 2) the Sierra) Nevada where 

D=LO 0 icneend =O. 40) Tires end Ss )ettne baci neend 

Range Province where b=o.4 km’and @a-I24 How,  Royecs al. 
(Mosby state that the silope associated’ with the hear 
tlow resulcs. Or theveasverm United Staves may nave 
Dreoadmanplicability im that Lb may be Considered. Che 
reference curve for normal continenval heavy, flow. | tne 
Basin and Range Province has an intercept of 1.4 9n-F2u.; 
O56 HebeUs @reater than the eastern Unived tater. 7 etnis 
Neat mMUst.come trom the) lower crust or upper mantles buy 
the lower crust. can De discounved as a soumce since (ais 
WOULG COMtlactewlun Seismic and @eccehnemical results. 
This leaves the Upper Manvle as, the source ol hel excess 
heat flux. Roy et al. (1968b) state that high tempera- 
(ees Ss een eane partial melting at shallow depths in. tine 


upper mantle would be consistent Wh ele 1eNiC wr eee Gi 
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and electrical conductivity data and that the TOC OF 
the partially molten zone is essentially an isotherm 
since the melting point gradient is only 2 to 4°C/kn. 

Even though the number of heat flow data loca- 
t2Ons is not as plentifulvas used) in other upper Manule 
investigation techniques in the continental United DLaALeS . 
enough Values exist in the western United States co 
support the concept of representing the surface or 
near suriace Of the Upper mantle a@e an s2sounerm that 
is undulating under the physiographic provinces in the 
wesverm Unived Staves. 

Warren et al. (1969) with heat flow results alone 
the same line as Schmucker (1969) were able to relate 
these heavy flow vVarlavions wich Schmucker's 
electromagnetic anomalies by showing that both sets of 
data were consistent with a 1200°C isotherm at a depth 
between 50 and 200 km under Southern Arizona and between 
150 and 400 km under western Texas. Warren's isotherm 
was determined by downward continuation of heat flow 


data in two dimensions Using Poisson's equation: 


Seay ee ee) 
2 2 K 
dX OZ 


where @ and A are the temperature and heat production 
at the point (x,z). He @ssumed the surface of the carth 


also to be an isotherm. The results of this process are 
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Pigure )5.. / Heat flow and conductive models 


(itex Warren ct alow 1960). 
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Pigure 5.38 Heat flow and geomagnetic depth-sounding 
and magnetotelluric stations (Warren et 
al. (960). The three broken lamess show 
Dae =e pOrOLiMa re POs tlOlem en sie cua ees 
and 4 ef the 1960/ array in the present 
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Figure 5.9 Model structure for results of 1967 
acray (Porathver seal, 1970).  Repetinon 
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Shown in figure 5.7 along with the results of Schmucker 
(1969). Estimates can be made of the isotherm which 
gives Dest Correlation of heat flew to magnetic varia- 
utons alone a line, “The resulting temperature estimates 
are in agreement with estimates from OUNer COnslderawsons 
within 300°C. 

Figure 5.8 reproduced from Warren et al.. (1969), 
shows the distribution of heat-Tlow values in 4 reeien 
which overlaps most of the area covered by our 1967 
erray. lanes «2, 3 and 4 of the array are andiicested. 
the correlavion with the conductive structures (figure 
bro is Close. The Rasim and Ranse Province ana coucnern 
ROocKwveSs are how, the?Texas Foreland and Colorado Plaveau 
eold. Only two heav-hiow Gdeverm nations sussesis es lower 
Upper Manvle vemperature under vhe Colorado: lavesau. 
where Warren et al. obviously assume a correlation of 
heave tlow witb viessuruc lures. Melee PNeClGD wer Ochce 
Oflassteo Struckure an@ the Conducaivity Under the 
boundary between, the Basin and Range and the western 
edge of the Colorado Plateau bears out and extends the 
heat-flow picture. The ridge on top of the step at the 
Wasatch Front (figure 5.9) is not shown in the heat-flow 
results, probably because the latter are TOO. ews 

An alternative presentation of the correlation 
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PaSpure 5.10. Here 2 provisional heat-flow contour map 
by Roy, Blackwell and Decker (personal communication) 
is superposed on the unseparated Z component Fourier 
spectral amplitude map at T=60 minutes for the sub— 
SLormerwel dor September 1, 1967. The magnetic fields 
indicate the local currents under the Southern Rockies 
and Wasatch Front, west of the maxima in Z. The 
tracking of the magnetic field with the heat-—flow 
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Seismic, heat-flow and geomagnetic depth sounding 
wesulvs: agree in Showine wnat the Basim and Rance sFro— 
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wemlperavures, ab a Given depth, Chan That under tne 
Groeau Plains (or tne Colorado Plaveau.. Masnevevelimiric 
resulLue oy owl et Glo?) alse: fan ai VoOmpmtor~ Lerurer 
Figure 5.9 presents the essential results of vine 1967 
GDS array study discussed in Chis tinesits.. Pie moder. 
phere eiven 1S, Of @OUrse, NOW Unique. “ip ts, however, 
consistent with all of the geophysical information, 
though the GDS array study resolves devaile nog revealed 
by the present seismic and heat-flow data. 

Other large GDS array studies have been made in 


1968 and 1969, to extend the investigation of upper 
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mantle conductivity structure southward to the Mexico-— 
U.S. border and northward to the Trans-Canada Hoshway. 
When results are available from these later arrays it 
SHOU De pessible v6 sev 1amits TOvtine structures 
relavine tnose here discusseu vo tie East Paci ic ase 


and to those under the Northern Rockies. 
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